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In  this  dissertation  a  formal  kinetic  theory  is  used  to  cast 

the  line  shape  function  into  a  form  that,  while  similar  to  the 

"unified"  theories  of  Smith,  Cooper,  and  Vidal  and  of  Voslamber, 

does  not  introduce  some  of  the  usual  approximations.  The  resulting 

line  shape  function  explicitly  includes  the  initial  correlations 

between  the  atom  and  perturbers,  and  also  demonstrates  the  natural 

separation  of  plasma  mean  field  and  collisional  effects.  The 

classical  path  and  no-quenching  approximations  are  discussed  and 

ultimately  employed;  however,  they  are  not  required  in  the  formal 

development.  The  weak  coupling  limit  is  considered  as  a  systematic 

approximation  to  the  formal  results.  It  is  shown  that  different 

ways  of  applying  this  limit  lead  to  different  expressions  for  the 

memory  operator,  some  of  which  correspond  to  existing  theories. 

One  approximation  is  considered  which  systematically  incorporates 

the  effects  of  electron  correlations  within  the  framework  of  a 

unified  theory.  In  addition,  a  practical  approximation  suitable 

for  a  strongly  interacting  plasma  is  discussed. 
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CHAPTER  I 
PLASMA  UNE  BORADENING 

I-A 


For  a  large  range  of  temperature  and  density,  most  of  the  radiation 
emitted  by  gaseous  plasmas  is  due  to  atomic  transitions.  If  individual 
radiators  are  considered  to  be  isolated  and  stationary,  then  the 
width  of  spectral  lines  will  only  be  due  to  natural  line  broadening. 
On  the  other  hand  if  the  radiators  move  and  also  interact  with  each 
other  and/or  the  plasma,  pressure  and  Doppler  line  broadening  must  be 
taken  into  account.  In  a  plasma  in  which  a  significant  percentage  of 
the  particles  are  ionized  the  dominant  factor  in  pressure  broadening 
will  be  Stark  boradening,  the  understanding  of  which  requires  the 
study  of  both  atomic  physics  and  the  many-body  physics  of  plasmas. 

Plasma  spectroscopy  has  been  shown  to  be  a  particularly 
useful  diagnostic  technique  for  laboratory  and  astrophysical  plasmas. 
This  is  due  to  the  fact  that  radiation  emitted  by  the  plasma  acts  as 
a  non«interfering  probe,  in  other  words,  it  is  not  necessary  to  disturb 
the  radiating  system  in  order  to  measure  the  spectrum  of  the  emitted 
radiation.  Impetus  exists,  therefore,  to  calculate  theoretically 
the  emitted  (absorbed)  line  shapes.  The  atomic  physics  of  light 
atoms  is  well  understood,  so  it  is  the  manyfcody  theory  that  presents 
the  chief  obstacle.  In  Chapters  I  and  II  of  this  dissertation  we 
discuss  some  of  the  recent  approaches  to  this  problem;  in  Chapters  III 
and  IV  an  alternative  approach  based  on  kinetic  theory  is  proposed 
and  developed. 


If  an  excited  atom  or  molecule  were  alone  in  a  radiation  field 
which  initially  contains  no  photons*  the  excited  electrons  could 
make  a  transition  to  a  state  of  lower  energy  thereby  emitting 
electromagnetic  radiation  with  a  frequency  approximately  equal  to 
the  energy  difference  between  the  initial  and  final  states*  in  units 
of  "ft.  In  order  to  obtain  a  power  spectrum  for  such  a  system*  we 
consider  an  ensemble  of  excited  atoms  or  molecules  and  take  an  ensecble 
average  of  the  radiation  emitted.  The  power  spectrum  thus  obtained 
in  composed  of  several  very  narrow  emission  lines  at  the  characteristic 

natural  frequencies  of  the  atom  or  molecule  in  question.  The  theory 

o 

describing  these  natural  lines  is  discussed  by  Heitler  who  points  out 

that  the  widths  of  optical  lines  for  atoms  are  on  the  order  of  lO""4 

Angstroms. 

3 

A  calculation  of  a  line  shape  that  includes  Doppler  broadening 

requires  a  knowledge  of  the  velocity  distribution  of  the  radiating 
atoms.  This  subject  along  with  a  very  thorough  discussion  of  the 
general  line  broadening  problem  is  given  by  Griem. 

When  these  excited  atoms  or  molecules  are  exposed  to  time  dependent 
perturbations  due  to  their  interaction  with  other  particles  of  the 

gas*  degeneracies  of  their  excited  state  energy  levels  may  be  removed 

4 

and  the  half-life  for  a  given  transition  may  be  appreciably  altered. 

This  is  called  pressure  broadening.     In  this  case  the  power  spectrum 
for  the  radiators  in  question  will  show  a  redistribution  of  frequencies. 
A  typical  power  spectrum  contains  broad  lines  which  are  usually  centered 


near  characteristic  natural  frequencies;  hence  one  says  that  the 
natural  lines  have  been  broadened  and  perhaps  slightly  shifted  by 
the  perturbations. 

If  significant  numbers  of  particles  are  ionized  then  the 
strongest  of  the  pressure  broadening  mechanisms  will  involve  the 
interaction  of  ions  and  electrons  with  the  radiating  particles. 
Since  electric  fields  are  involved  this  interaction  is  usually 
called  Stark  broadening.    For  atoms  in  many  laboratory  plasmas 
this  effect  is  quite  important,  typical  broadened  widths  being  on 
the  order  of  10 3  times  larger  than  the  natural  width.     In  fact, 
the  broadening  of  spectral  lines  in  fairly  dense  plasmas  (e.g.  10 16 
per  cm3,  104  degrees  K)  is  due  almost  entirely  to  the  Stark  effect 
which  is  also  several  orders  of  magnitude  larger  than  the  broadening 
due  to  neutral  atom  pressure  broadening  mechanisms. 

For  typical  plasma  conditions  dealt  with  in  this  dissertation 

3 
Doppler  boradening    will  only  be  significant  near  the  very  line 

center.    However,  this  effect  is  easily  included  in  the  formalism 

developed  in  Chapters  III  and  IV. 

I-C    S.tgrK,  proadeniag.  JP.  PJtefflS, 
There  are  two  boradening  agents  responsible  for  the  Stark  effect 
in  ionized  gases:     ions  and  electrons.    The  broadening  caused  by  each 
of  these  is  considerably  different  due  to  the  difference  in  their 
velocity  distributions.       In  order  to  illustrate  this  difference, 
we  note  that  the  length  of  time  which  is  of  importance  in  line  broadening 
is  the  half-life  for  the  excited  state  of  the  atom.     In  a  few  half- 
lives  we  may  assume  that  an  atom,  originally  in  an  excited  state. 


has  radiated,  hence  any  process  that  takes  many  half-lives  will  be 
almost  static  from  the  point  of  view  of  the  excited  atom.     It  is 
just  this  fact  that  enables  us  to  distinguish  electrons  from  ions. 
The  ions,  being  much  heavier  than  the  electrons,  move  more  slowly 
and,  for  most  plasma  problems,  the  distribution  of  the  ions  does 
not  change  appreciably  during  a  few  half-lives.    This  is  the  basis 

for  the  "quasistatic"  or  "statistical"  approximation  which  was 

5  6 
developed  by  Holtsmark.  *      The  meihod  of  treating  the  electron 

collisions  represents  an  example  of  the  opposite  extreme,  in  which 

a  process  takes  a  very  small  fraction  of  a  half-life  and  may  be 

regarded  as  instantaneous.     An  approximation  based  on  this  limit, 

1  7 
the  impact  approximation,  *      has  been  used  to  treat  the  electrons 

since  their  great  speed  causes  most  of  them  to  pass  rapidly  by  the 

atom  producing  collisions  of  very  short  duration  compared  to  the 

excited  half-life.    These  two  limits  will  be  discussed  further  in 

Section  I-E  of  this  dissertation. 

i-d  Tte.&pp,  Shape 

The  radiation  spectrum  of  a  quantum  system  is  determined 

experimentally  by  nfiasureroent  of  the  power  radiated  per  unit  time 

per  unit  frequency  interval,  averaged  over  the  polarization  and  the 

direction  of  radiation.  Since  this  quantity  has  been  derived 

many  times  in  a  variety  of  ways,  its  derivation  will  not  be 

included  here.     The  power  radiated  when  a  particle  makes  a  spontaneous 

dipole  transition  from  one  quantum  state  to  another  is 

4  a)2 
P(u))=37§-KJ,  (I-D-1) 


where  I(oj)  is  usually  referred  to  as  the  line  shape.  It  is  often 
convenient  to  write  I  (to)  in  terms  of  its  Fourier  transform  $(t), 

1(0))=  -J-  f  dteiaVYt*(t)  (I-D-2) 

^  '  mm  00 

«,(t)=rr{d.T(t)(Pd)T+(t)}.  (I-D-3) 

The  damping  factor,  y  ,  v/hich  represents  the  effect  of  the  natural 
half  width  of  the  spectral  line  of  interest  was  introduced  in  an 
ad  hoc  manner  by  Smith.      While  the  effect  of  this  term  is  negligible 
over  most  of  the  spectral  line,  it  is  a  convenient  quantity  to  introduce 
since  it  will  insure  the  convergence  of  certain  indeterminate  integrals. 
Tr{•••}  denotes  the  trace  operation  taken  over  the  plasma-radiator 
system,     d  is  the  dipole  operator  for  the  radiator  and  T(t)  is  the 
time  development  operator  for  the  system  .     The  latter  satisfies  the 
usual  equation  of  motion; 

i-hft  T(t)=OT(t)  (I-D-4) 

and  can  be  written  formally  as 

T(t)=e       *     ;  T(0)=lj  (i-D-S) 

where  H  is  the  Hamiltonian  operator  for  the  system,  the  form  of  which 
will  be  specified  later.  The  density  matrix  for  the  plasma  radiator 
system  when  it  is  in  thermal  equilibrium  is  given  by 

P-e   /Tr  ie       i 

3=1/^  (l-D-6) 

where  IC  is  Boltzmann's  constant  and  T  is  the  temperature. 


The  Fourier  transform  of  the  line  shape.,  $(t)M   can  be  interpreted 
as  the  autocorrelation  function  for  the  amplitude  of  the  wave  train 
of  the  radiation  emitted  when  a  radiator  makes  a  spontaneous  dipole 
transition.   An  important  property  of  $(t)  is  that 

*(-t)=I*(t)l*  (I"°"7) 

which  allows  us  to  write  I(cj)  as  a  Laplace,  rather  than  a  Fourier 
transform, 

l(co)=i-R  fdteiaVYt*(t).  <I'*-8) 

■n    eJQ 

This  is  an  important  difference  since  the  initial  value  of  any  equation 

of  motion  may  be  specified  more  easily  at  t=0  rather  than  in  the 

difficult  limit  t->--°°. 

The  Hamiltonian  operator  H  is  taken  as  the  sum  of  three  terms, 

H=H  +H  +HT  (l-D-9) 

a     p    I 

where  H  is  the  Hamiltonian  for  the  isolated  radiator,  H  is  the 
a  P 

Hamiltonian  for  the  isolated  plasma,  and  Hj  is  the  Hamiltonian  for 

the  interactions  between  the  radiator  and  the  particles  of  the 

plasma.  H  and  HT  can  be  further  subdivided  into  the  contributions 
P     I 

from  electrons  and  ions 


H  =a  +H.+H. 
p    e     1     le 

H=tt     +H.      .  (I-D-10) 

I     ea    1a 

4 

It  is  convenient  now  to  define  the  Liouville  operator,  L,  given  by 

Lf=    £tH,f]j      L^tH^f]}   etc.  (1-43-11) 


We  can  now  define  the  Liouville  time  development  operator, 

-iHt    iHt 
e-iLtf=e  f   fe"n  -r^jfxt^)  u  (I-D-12) 


With  this  definition  the  autocorrelation  function,  given  by  I-D-3, 
becomes 

$(t)=Tx{^«e"iLtp^}.  (I-D-13) 

Equation  I-D-13  is  now  seen  to  satisfy  the  Liouville  equation,  a 
property  that  will  be  utilized  in  Chapter  III. 

I-E  Model  For  the  Plasma 

In  order  to  evaluate  equation  I-D-2,  additional  approximations 
are  required.  There  have  been  a  number  of  different  theoretical 
approaches  to  the  problem  of  line  broadening  in  plasmas,  each  of 
which  has  employed  a  different  set  of  assumptions.  In  this  section 
we  will  discuss  a  few  of  the  approximations  that  are  common  to  most 
of  the  theoretical  approaches.  We  will  further  consider  what 
restrictions  these  place  on  the  model  of  the  plasma. 

No  Quenching  Approximation 

One  important  approximation  invoked  by  most  theoretical 

calculations  of  a  line  shape  is  the  so-called  no-quenching  approxi- 
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mation.  This  has  been  discussed  in  detail  by  many  authors  *  *  and 

we  will  merely  outline  its  implications  here. 

Any  complete  set  of  states  may  be  used  to  evaluate  the  trace 

operation  and  the  matrix  elements  contained  in  equation  I-D-2. 

However,  it  is  usually  most  convenient  to  use  the  complete  set 


formed  by  the  direct  product  of  the  free  radiator  states  denoted 
by  | a>  and  the  free  plasma  states  denoted  by  | a>* 

|  aa>=|  a>  |  a>  •  (WML) 

These  states  satisfy  the  following  eigenvalue  equation 

Ha|a>=£a|a> 
Hp|a>^Ja> 

(H  +H  )|aa>=(E  ■•£   )|aa>  (l-E-2) 

a    p  '  a    a   ' 

where  E  and  E  are  the  energy  eigenvalues  for  the  free  radiator 
a     a 

and  free  plasma  respectively.  The  initial  and  final  states  of  the 
radiator  will  be  denoted  by  |i>  and  | f >  respectively. 

The  autocorrelation  can  be  rewritten  by  expressing  the  trace 
as  a  sum  over  products  of  matrix  elements, 

$(t)=    z      <aa[ar|a-a'><a'a'|T(t)!a"a"><a"a"|p3'|a"'ci"->x 
a, a  ,  a  ,  a 

x<a— a"-|T+(t)|aa>.  (I-E-3) 

We  shall  be  interested  in  the  radiation  that  results  when  an  atom 
spontaneously  decays  from  some  excited  state  of  principal  quantum 
number  n  to  some  state  of  lower  principal  quantum  number  n'.     For 
low  lying  levels  of  hydrogenic  atoms  it  is  commonly  assumed  that 
radiationless  transitions  are  improbablej  hence,  T(t)  and  Tf(t) 
will  be  assumed  to  have  matrix  elements  only  between  states 
of  the  same  principal  quantum  number.    This  restriction,  which 
is  the  no-quenching  approximation,  allows  us  to  rewrite  I-E-«3 


*(t)=£<fa|2|  ia><ia|T(t)|  iV><i*a*|p3|  f* a"><ra"|T+(t)|  fa>         (!•*■*) 

where  we  have  used  the  fact  that  cl,  a  pure  atomic  operator,  will  be 
diagonal  between  free  plasma  states* 

For  simple  radiating  systems  for  which  the  atomic  physics  is 
well  understood  the  real  problem  that  must  be  considered -is  the 
statistical  average  over  perturber  states.     In  general  this  is  an 

extremely  difficult  problem,  but  one  v/hich  can  be  simplified  by 

1  7 

employing  the  classical  path  approximation.     In  this  approximation 

the  wave  packets  of  the  perturbers  are  assumed  to  be  small  enough 
so  that  they  do  not  overlap  either  with  each  other  or  with  the 
radiator.  This  allows  us  to  view  the  perturbers  as  point  particles 
traveling  in  classical  trajectories  and  interacting  with  classical 
potentials.  The  effect  of  this  approximation  is  to  replace  p, 
defined  by  equation  I-D-6,  by  its  classical  analog  and  correspondingly 
to  replace  the  trace  over  perturber  states  by  integrals  over 
perturber  coordinates.  These  simplifications  are  probably  valid 
for  ions  except  at  extreme  densities  and  low  temperatures,  but 
since  electrons  are  much  lighter  and  their  wave  packets  much  larger 
they  could  present  a  problem.  Several  authors  have  shown,  however, 

that  over  much  of  the  temperature  and  density  ranges  considered  by 

1  7 

experiment  the  electrons  will  indeed  behave  like  classical  particles.  * 

Factoring  The  Initial  Density  Matrix 

An  additional  approximation  that  is  made  throughout  much  of 


JLU 

the  literature  *  *  *     is  the  factorization  of  the  density  matrix, 

P=Papipe*  (l"C«6) 

where  p  ,  p.   and  p    represent  the  density  matrices  for  the  radiator*  the 
ions,  and  the  electrons  respectively.    This  factorization  implies  that 
the  interactions  among  the  ions,  the  electrons,  and  the  radiator  have 
been  neglected.    Static  correlations  between  the  ions  and  the  electrons 
are  not  ruled  out  entirely,  however,  because  they  can  be  partially 
accounted  for  by  replacing  the  ion-ion  interactions  that  appear  in  p 
by  some  effective  interactions  which  attempt  to  account  for  electron 
screening  of  the  ions.     On  the  other  hand,  the  perturber-atom  interactions 
are  neglected  in  p;  the  effect  of  this  exclusion  will  be  discussed 
later.     Using  this  factored  form  of  the  density  matrix,  the  average 
over  perturber  coordinates  in  equation  1-0-13, 

$(t)^rr{3-e"iLtp3},  (I-E-6) 

can  be  divided  into  a  trace  over  atomic  coordinates  and  a  trace  over 
perturber  coordinates 

$(t)=Tra{d-Trie{e"lLtp.pe}pad}.  (l-C-7) 

7 

Bar anger  has  shown  that  this  approximation  is  equivalent  to  the 

neglect  of  back  reaction;  that  is  to  say  the  trajectories  of  the 
perturbers  are  assumed  to  be  unaffected  by  their  interactions  with 
the  radiator. 

Until  recently  most  theories  of  line  broadening  have  assumed 
that  the  ions  are  static:  they  are  assumed  to  be  so  heavy  that  their 


distribution  is  not  appreciably  altered  during  the  radiative 

lifetime  of  the  atom.  This  is  equivalent  to  assuming  that  the 

4 
ions  are  infinitely  massive.  Smith  has  shown  that  this  approximation 

implies  that  the  kinetic  energy  part  of  the  ion  Hamiltonian  will 

commute  with  the  potential  energy  part.  The  effect  of  this  commutation 

is  that  the  free  ion  Hamiltonian  does  not  appear  in  the  Liouville 

time  development  operator  in  equation  I-£-7.  This  allows  us  to 

write  equation  I-D-8  in  the  following  form 


1(0))  = 

where 


deP(e)j(o)e)  (l«£-8) 


j(coe)=  *-R 


dteia)tTr  {d'e-iL(e)tp  d}  (l-E-9) 


ae         ae 


and  where  P(0  is  the  static  microfield  distribution  calculated  by 
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Hooper.     It  will  be  noted  that  all  of  the  dependence  on  the 

ion  field  variable  is  contained  in  the  Liouville  operator,  L(e). 

The  associated  Hamiltonian  is  then  given  by 


H(e)=H  +ee-R+H  +H  (l-E-10) 

a     e  ae 

where  the  ion  radiator  interaction  is  given  in  the  dipole 

approximation.    R  is  the  position  vector  of  the  atomic  electron 

and  e  is  the  electric  field  strength  at  the  radiator.     For 

convenience  in  notation  we  will  define 

H  (e)=H  +et-$ 
a  a 

which  allows  us  to  write 

H(e)=H  (e)+H  +H  (i-E-ll) 

3      3©   6 
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L(e)=La(0+Lae+Le.  (I-&42) 

As  a  result  of  the  static  ion  approximation  we  are  able  to 
concentrate  on  the  electron  broadening  of  an  atom  placed  in  an 
external  field;  that  is,  we  determine  J(<de). 

The  effect  of  ion  motion  has  been  shown  to  be  important  only 

4 
in  the  very  center  of  a  spectral  line,  but  there  are,  nevertheless, 

certain  physical  situations  where  it  is  observable.  Several  authors 
have  discussed  this  problem*   "  '  It  will  be  ignored,  however,  through- 
out the  remainder  of  this  dissertation,  and  the  ions  will  henceforth 
be  assumed  to  be  static. 


CHAPTER  II 
UNIFIED  THEORIES 

II-A     Introduction 


Before  the  late  1950*s  it  was  thought  that  all  Stark  broadening 
resulted  from  ions  since  the  fast  moving  electrons  were  expected  to 

have  no  net  effect.  '      It  was  first  shown  by  Kolb  and  Griem      in 

17 

1958  and  by  Baranger        that  electron  dynamics  would  be  important 

and  result  in  a  considerable  amount  of  broadening.    Exact 
calculations  of  line  profiles  that  include  the  effect  of  the  electrons 
are  not  in  general  possible*  hence  a  number  of  theoretical  approaches 
have  appeared,  each  of  which  uses  a  different  set  of  approximations. 

All  line  broadening  theories  may  be  divided  into  two  broad 
categories,  those  which  are  fully  quantum  mechanical  are  those 
which  make  the  classical  path  assumption  for  the  perturbers.     Several 
fully  quantum  mechanical  theories  have  been  developed*  however  the 
only  successful  calculations  of  entire  line  profiles  have  been 
based  on  theories  in  which  the  classical  path  approximation  has  been 
made.     In  the  remainder  of  this  chapter  we  will  discuss  only  these 
theories.     A  fully  quantum  mechanical  treatment,  valid  over  the 
entire  line  shape,  would  be  preferable,  of  course,  but  this  has  not 
yet  been  possible.     It  should  be  noted  that  the  formalism  developed 
in  Chapter  III  of  this  dissertation  is  not  limited  to  the  classical 
path  approximation  and  might  be  used  as  a  starting  point  for  an  all 
order,  fully  quantum  mechanical  calculation. 
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The  first  line  shape  calculations  that  included  electron 
broadening  realistically  were  the  so-called  impact  theories  of  Kolb 
and  Griem,  and  of  Baranger.  While  these  approaches  employ  a  factored 
density  matrix  and  make  the  static  ion  and  classical  path  approximations, 
they  also  make  the  impact  approximation  and  the  completed  collision 
assumption.  The  completed  collision  assumption  treats  the  electrons 
as  moving  so  fast  that  collisions  can  be  considered  instanteous. 
Baranger  has  shown  that  this  is  equivalent  to  a  Markov  approximation 
that  leads  to  a  considerable  simplification.  Also  included  in  the 
impact  theory  is  the  impact  approximation.  This  assumes  that 
close  collisions  between  atom  and  electrons  will  occur  one  at  a 
time  and,  hence,  furiher  simplifies  the  calculation  by  considering 
only  binary  electron  atom  collisions. 

The  early  impact  theory  calculations  also  took  the  collision 
operator  (Chapter  III)  to  second  order  in  the  electron-atom 
interaction.  Any  theory  that  makes  this  approximation  is  called 
a  second  order  theory.  It  should  be  pointed  out  that  not  all 
second  order  theories  are  impact  theories:  Smith  and  Hooper  in 
1967     developed  a  fully  quantum  mechanical  second  order  theory 
without  the  completed  collision  or  impact  approximations. 

It  will  be  observed  in  Chapter  IV  that  second  order  theories 

21 

are  valid  only  in  the  line  center,  while  one-electron  theories, 

which  result  from  expanding  J(we)  in  powers  of  the  inverse  frequency 
separation  are  valid  only  in  the  line  wings.    Recently  several  unified 
theories  have  been  developed  that  are  valid  both  in  the  line  wings 
and  in  the  line  center. 
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II-6     adfifij  Theory  of  Vidal.  Cooner  and  Sqjth 

The  first  successful  unified  theory  calculations  were  carried 
out  by  Vidal,  Cooper,  and  Smith  (VCS).22*23    Although  they  used  the 
classical  path  and  static  ion  approximations  together  with  a 
factored  density  matrix  they  included  all  orders  of  the  electron 
atom  perturbation. 

The  starting  point  for  the  VCS  calculation  is  the  electron 
broadening  line  shape,  given  by 

JU)=  K  r<*teiuV6^L(e)iW}  (II"B-1) 

tt    ejp  ae  d  e 

where  the  approximation  that  Pae=PaPe  haS  besn  utilized.    The  form 
of  II-B-1  can  be  simplified  by  separating  the  trace  over  atomic 
coordinates  from  that  over  electron  coordinates  and  by  performing 
the  Laplace  transform: 

24 

VCS  then  employ  a  projection  operator  technique  developed  by  Zwanzig 

or  96 
and  by  Mori    *       to  obtain  the  result: 

J(ue)=  LR  Tr  {         *TT   '   pa>  (II-B-3) 

11    e     a  anLa-M(w)     3 


where 

M(a))=-i 


°°dteiAa)0Pt<L     (t)G(t)L     > 
0  ea  ea 


G^Kxp{-  Jr 


rt 


dt'(l«P)L,Jt')} 
(e). 


0  ea 


A>^La(£);     L^^^^K/^^K  d«-4) 
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T  is  the  time  ordering  operator  and  P  is  the  projection  operator 
defined  by 

Pf=PeTre(f}.  (IIH3-5) 

In  equation  II-6-3  all  of  the  complicated  dependence  on  the  N  electrons 
has  been  transferred  into  the  effective  atomic  operator,  M(o)).  This 
is  a  particularly  useful  functional  form  for  the  line  broadening 
problem  since  M(w)  can  be  interpreted  as  a  frequency  dependent  width 

g 

and  shift  operator. 

Equation  II-*-4  is  still  not  a  form  which  is  amenable  to  calcu- 
lation.    One  means  of  approximating  this  equation  is  to  set  the  time 
development  operator  in  II-6-6  equal  to  1.     This  leads  to  the  so- 
called  second  order  theories  which  receive  their  name  from  the 
fact  that  M(w)  is  second  order  in  the  electron-atom  interaction. 
Many  line  broadening  calculations  have  been  made  using  this 
approximation:     from  the  early  calculations  of  the  impact  theory 
to  a  fully  quantum  mechanical  treatment  by  Smith  and  Hooper.  ^^  ' 
However,  this  approximation  breaks  down  in  the  line  wings.     In  other 
words  when  L      and  therefore, M(w),  is  large,  it  is  incorrect  to  stop 
at  finite  order  in  an  electron-atom  perturbation  expansion. 

Since  the  unified  theory  is  supposed  to  be  valid  in  the  line 
wings  as  well  as  in  the  line  center,  setting  the  time  development 
operator  equal  to  1  is  not  allowable.     VCS  showed  that  to  lowest 
order  in  the  density  the  time  development  operator  could  be 
approximated  in  the  following  way; 

ft 

G(t)=T       {-£      «'•""»*(  WOsLfjfc11**} 

exp      Ti  j  q  j 

*nT      UJLPdt-L  a(j.t-)}  (ll-B-6) 

j  exp     -fi  J0        ea 


with  M(w)  given  by 

M(a))!=^NfdteiMbPt<L     (l,t)T       {-> 
J  ea  exp      Ti 


dtT.    (l,t)}L    (l).     (II-B-7) 


The  average  is  performed  over  the  coordinates  of  particle  1  only. 
It  should  be  noted  that  I 1-6-7  is  all  orders  in  the  electron-atom 
interaction.     It  was  shown  that  in  this  approximation 

(l-P)L     (l)-L     (1).  (II-B-8) 

Co  6  a 

One  important  limitation  imposed  by  taking  only  the  lowest  order 
term  in  the  density  expansion  is  that  electron-electron  interactions 
are  omitted  both  from  p  and  L    (j,t).    VCS  partially  correct  for 

Q  3 

this  by  assuming  that  the  electrons  can  be  replaced  by  shielded 
quasi-particles,  in  other  words  they  assume  that  the  electric 
field  in   e  «R  is  shielded.     We  will  see  in  the  next  section 
how  Capes  and  Voslamber  included  electron-electron  correlations 
in  their  theory. 

Another  VCS  approximation,  the  neglect  of  time  ordering  from 
equation  II-6-7,  has  been  discussed  in  detail  by  Smith,  Cooper,  and 
Roszman,       and  it  will  not  be  further  considered  here.     Integrating 
by  parts,  equation  II-C-7  becomes 

M(w)~ihAu)  Jdte^V^expI-j- [  dt'L     (l,t')}-l>Aa>     .  (II-B-9) 

opj  n  j  q         ea  op 

Finally  VCS  assume  that  the  time  dependence  of  R  can  be  neglected 

in  the  expression  for  L  (l,t);  hence 

e  a 

L  (l,t)=eS-e"iI^(l)tt(l)ei^(l)t.  (II-B-IO) 

6  3 
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Instead  of  actually  performing  a  spatial  integration  over  the 
time  dependent  shielded  electric  field,  •'4MX'te(l)e11^   * 
VCS  replace  it  by  an  unshielded  one  and  cut  off  the  resulting 
integral  over  the  free  particle  trajectories  at  the  Debye  sphere. 
It  is  interesting  to  note  that  if  the  upper  limit  of  the  integral 
in  the  exponential  of  II-B-^  is  extended  to  infinity,  the  time  develop- 
ment operator  in  the  interaction  representation  goes  over  to  the 
S-matrix  and  we  regain  an  all  order  impact  theory.  This  is  entirely 
equivalent  to  the  completed  collision  assumption.  Thus  the 
expression  for  M(u)  finally  evaluated  by  VCS  was 

,t 

c 
0 


t 
M(u)=^4nAu)0  [dte^oP^expC-M  dt'Lea(r4  t')}-l>Aa>op.   (H-B-ll) 


The  resulting  matrix  was  inverted  to  yield  an  expression  for  the 
line  shape. 

The  numerical  results  obtained  by  this  formulation  of  the 
line  broadening  problem  has  been  shown  to  agree  well  with  most 
currently  available  experimental  results.  Nevertheless,  the 
approximations  that  have  gone  into  this  development  are  not  entirely 
transparent.  The  density  expansion  was  truncated  in  a  way  that 
precludes  the  possibility  of  going  to  higher  order.  In  addition, 
there  is  no  way  of  examining  the  validity  of  the  way  in  which 
electron  correlations  were  included. 

It  is  also  possible  to  formulate  a  systematic  kinetic  theory 
approach  to  line  broadening.  This  method,  as  developed  by  Capes 
and  Voslamber,  *   will  be  discussed  in  the  next  section. 
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II-C    Theory  of  Capes  and  Voslamber 

A  systematic  treatment  of  the  line  broadening  problem  may 
be  effected  by  viewing  the  radiating  atom  and  the  N  electrons  as 
an  N+l  particle  system  obeying  a  Liouville  equation.     Capes  and 

97  OQ 

Voslamber"  *   have  developed  such  a  theory  that  included  electron 
correlations  in  a  more  systematic  fashion.  This  section  will  outline 
their  approach,  indicating  its  advantages  and  where  it  differs  from 
the  VCS  theory. 

Capes  and  Voslamber  make  several  of  the  usual  approximations 
including  the  classical  path  approximation,  the  neglect  of  perturber- 
atom  interactions  from  the  density  matrix,  and  the  static  ion 
approximation.  They  do  not,  however,  neglect  electron  correlations, 
and  their  theory  offers  an  understanding  of  what  effect  these  correlations 
will  have  on  the  line  shape. 

As  in  the  VCS  theory,  Capes  and  Voslamber  start  with  the 

equation 

ft 
ir  ~  e 


J~U)=  LR    fdte^r  {1-Tr  {•*iL(e)tp >  3}  (Il-Ol) 


a  el  re^a 


which  they  reduce  to  a  simpler  form  by  utilizing  the  fact  that  the 
integrand  obeys  a  Liouville  equation: 

(|r-+iL  +i£L     (jlK^l'-NjtHlVV^'  •2--)D(a,l'"N;t)=0 
8t        a     j  ea  .•  J     j     m    ■+         -► 


3xj       8V  j 


Tk1 
where  we  have  defined 


6(a,l— N|t)«VNe"iL(E)tP  d  3,  (II-OO) 


e^a 
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and  where  we  have  explictly  made  the  classical  path  approximation. 
Taking  partial  traces  of  H-O-2,  a  3BGKY  hierarchy  of  kinetic  equations 
results 

(L  +iL  +i  Z  L     (j))&(a,l—«it)+  I  (v.-V  -  J.      I      -Ji  .  V>D(a,l—«|t) 
V3t         a     j=l  ea  j=l     J     J     mi=s+l   3x.  3v. 


—en 


where 


d(«+l)LM(««l)ff(a.l—«f«+l!t)-  j£-   d(s+l)  s    — s^-'^-D(a,l—sjt) 

i=l      9Xi       av.  (II-c-4) 


ea 


m 


ff(a,l— «it)arr6+x#..N{?(a,l—N3t)}. 

The  first  two  equations  of  this  hierarchy  are 


(II-C-5) 


(•3-  +iLa)U(ait)=-inJd(2)lf(a,ljt)- 

fr  +iL  +iL ■     (1))-HT  -V  )U"(a,lst)=-infd(2)L     (2)fj(a,l,2;t) 

oi         «      69  l      i  ea 


(II-06) 


-i 


EL 


m 


3<f), 


d(2)  -mfL*  ,a—  tf(a,l*2jt). 

■  3X  3V 

Since  J(we)  could  be  written  as 


(n«c-7) 


J(a)£)=  ^R    fdteiaStTr   {d*-tf(a;t) }, 


tt  e 


(II-C-8) 


Capes  and  Voslamber  then  solved  for  D(ajt).  To  do  so  they  assumed 
that  $"(a*l>2jt)  could  be  expressed  as  an  approximate  functional  of 
tf(ajt)  and  l?(a,lst).  This  resulted  in  a  pair  of  coupled  differential 
equations  which  were  solved  simultaneously.  Their  closure 
relationship  based  on  the  weak  coupling  approximation  was 


J?(a,l,2jt)=^a,l;t)Vp(2)+rT(a,l5t)VD(l)-^(a;t)v2(p(l2)-p(l)p(2))  (lI-C-9) 


<-± 


30 
where  the  reduced  density  matrix  is  given  by 


(l*»»s)=  ••• 


d(s+l)-«.d(N)o  . 


If  the  closure  relationship,  II-C-9,  is  substituted  into  the 
second  equation  of  the  hierarchy,  II--C-7,  and  if  only  those  terms 
which  are  lowest  order  in  an  expansion  in  the  electron -electron 
coupling  strength  are  kept,  then  II-07  becomes 

a  *     ->      ->  n    3MV|)f  9((>?1 

(|7+iLa+iLea(lK1^1)l5(a.ljt)».  J-V-    d<2)  ^■•(•-llt) 

d  V, 


1 


-inf0(Vl) 


d(2)Lea(2)D(a,2jt)-«f0(v1)Jd(2)Lea(2)g(12)p(2)S(a;t)s 

p  (12)=P  (Dp  (2)g(12)j     fQ  {vl  )^p  (l),  (n-c-10) 


where  the  symmetry  properties  of  the  interactions  have  been  used 
to  eliminate  some  of  the  resulting  integrals.  The  first  term  on 
the  right  hand  side  of  equation  II-C-10  is  recognized  as  the  Vlasov 
operator,  V(l).    If  the  Laplace  transforms  of  II-C-6  and  II-C-10 
are  solved  simultaneously,  we  find  that 


D(a,lsw)= _i Lea(l)D(aj  co)j      (n-C-ll) 

co-L -L  (l)+iv ,-V  -V(l) 
a  ea     1  1 

V(l)  is  again  the  Vlasov  operator,  and  L      is  the  electron-atom 
interaction,  statically  shielded  by  the  electron-electron  pair 
correlation  function.     If  equation  II-C-11  is  substituted  into 
the  right  hand  side  of  the  Laplace  transform  of  II-05  then  we  get 


SUi«">°  „-L  irit.)  >*  »««») 
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where  the  memory  operator,  M(o>),  is  given  by 
M(u  )=~in 


.s 


d(1)L-(1)  l  i  (i)  r ; ,(,) tJ1)'      uwH») 

a    ea  i     j 

Using  a  technique  similar  to  the  one  presented  in  Appendix  F  Capes 
and  Voslamber  showed  that  this  equation  reduces  to 

M(u))=-iN[d(l)  J-  fdteiAwoptLD  (1/t)u(t)L*  (i) 
J  2tt  J  ea  ea 

where  L ■     (l*t)  is  a  dynamically  shielded  interaction  and  U(t)  is 
ea 

the  interaction  representation  time  development  operator.    This 
result  is  identical  to  that  of  VCS  except  for  the  shielding  which 
appears  in  the  electron-atom  interactions.    Note,  however,  that 
the  electron-electron  interaction  appearing  in  U(t)  is  not  shielded 
while  the  two  interactions  appearing  around  it  are.     Thus  it  is  not 
possible  to  further  simplify  II-C-14  by  performing  an  integration 
by  parts  as  was  done  in  equation  II-B-10. 

The  main  strength  of  the  approach  of  Capes  and  Voslamber  is 
the  fact  that  their  theory  includes  the  effect  of  electron 
correlations  in  a  way  which  is  preferable  to  the  ad  hoc  cut-off 
prdcedure  employed  by  VCS.     Its  main  weakness,  however,  is  the 
closure  hypothesis,  equation  II-07.     Capes  and  Voslanber  show 
that  this  is  closely  related  to  the  impact  approximation  (binary 
collision  approximation  for  electron-atom  interactions,  and  like 
the  impact  approximation,  there  is  no  clear  cut  way  of  improving 
upon  it.     In  the  next  chapters  we  will  close  the  hierarchy  in  a 
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way  in  which  approximations  may  be  more  systematically  made,  and 
we  will  suggest  an  approximation  procedure  which  leads  to  a  more 
inclusive  result  than  found  in  equation  II-C-14. 


CHAPTER  III 
FORMAL  KINETIC  THEORY  APPLIED  TO  LINE  BROADENING 

III-A    The  Hierarchy 

In  this  section  we  will  develop  a  kinetic  theory  of  line 
broadening  in  plasmas  similar  to  that  derived  by  Voslamber.     * 
The  formalism  will  be  developed  for  an  atom,  perturbed  by  an 
external  field,  and  immersed  in  a  one  component  plasma;  therefore 
the  static  ion  approximation  is  implicit.     However,  two  other 
frequently  employed  approximations  will  be  avoided:     the  classical 
path  approximation  and  the  neglect  of  electron-atom  interactions 
in  the  initial  density  matrix. 

From  equation  I-E-9  the  line  shape  function  for  electron 
broadening  is  given  by 

ir     e     ae        u)-L(e)  HaN      *  IHI-A-IJ 

This  can  be  written  in  the  alternative  form: 

~  t°° 

J(a,e)=*'ReJ  dteiwtTrae{d.  Jj-D(a,l—N,t) } 

JU)=  7HeTrae{3.  J-   D(a,l— Hj«)>  (III-A-2) 

where  we  have  defined  the  operator  D(a,l«»»Njt)  and  its  Laplace 
transform,  D( a, l#,,Nj(ij),  as  in  equation  H-O-3; 
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•zo 


D(U...Nlt)^e"iL(£)VJ 


D(a,l— Hjoi)^  ^j^y  p^cf  (HI-*-3) 

where  V  is  the  system  volume  and  p„     is  given  by 

paN*"e  A*™  }* 

H(e)=«a(e)+ZHe(j)+  E  H^+ZH     (j). 
5  i<j        j 

In  the  above  equation  H  (j)  represents  the  kinetic  energy  of 

i  J 
particle  j,  H  **  represents  the  interaction  between  particles  i  and 
ee 

j,  and  H     (j)  represents  the  interaction  of  the  atom  with  electron  j. 
ea 

The  Liouville  operator  corresponding  to  H(e)  is  given  by 

IFL  +EL  (j)+  E   L*f+     EL     (j)  (III-A-4) 

3  j  6         i<j  6e       3  ea 

where  here  and  henceforth  the  functional  dependence  of  L  md  L    on 

a 

the  ion  microfield,   e  ,  will  be  suppressed.    We  also  define  the  reduced 
functions. 


««,!...  6*t)=VSTr  +1..rf,{«"4Ltp.J}. 


s+l»»«N 


D(a,l-.Sj   )=VSTr  {  -i- p  J} 

s+l»»»N     oj-L  Ka 

p(a,l...s)=Trs+1{PaN>.  (III-A-6) 

Hence  the  expression  for  J(We)  can  be  written  as 

j(u)e)=^ReJ  dteia)tTra{d.D(ajt)}, 

j(coe)=  J-R  Tx  (d.D(aico)},  (III-A-<S) 
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and  thus  the  problem  is  to  determine  ft(a;t)  or  DCa;^). 

We  start  by  noting  that  the  operator  6(a,l««»Njt)  obeys  the 

31 

Iiouville  equation, 

(  §-  +iL)D(a.l...N;t)=0,  (III-A-7) 

from  which  we  can  now  generate  a  hierarchy  of  kinetic  equations 
by  taking  partial  traces; 

(  ~+iL(s))6(a,l"-s;t)~in  E  Tr  ..   „{(L  (c+D+L1*8*1) 
dt  ._,  s+l»»»N  N  ea      ee 

C(a,l...s+l;t)         (III-A-8) 

s      s 
LS=L  +  E  I   (j)+  E  L  (j)+  E   Lij 

a  j=l  e    j=i  ^    .<jiS  ee  (IIIWH?) 

where  the  thermodynamic  limit  has  been  assumed.  The  first  three 
equations  of  this  hierarchy  are  given  explicitly: 

(—•+iLa)0(a;t)=^-inTr1{Lea(l)&(a,lst)}, 

(§£  +iL(l))0(a,ljt)=-int2{(Lee(2>t|^)]5(a,l,2jt)}  , 

(ft  +iL(2))^(a,l,2jt)=-inl-3{(Lea(3)+Le3^Le3^0(a,l,2,3;t)}.     (m-A-lO) 

The  approach  followed  by  Capes  and  Voslamber  to  solve  for  D(a;t)  was 
to  use  a  closure  relationship,  expressing  D(a,l,2;t)  as  an 
approximate  functional  of  D(ajt)  and  D(a,l;t).  While  this  procedure 
led  to  a  closed  set  of  coupled  differential  equations  that  could  be 
readily  solved,  its  weakness  lay  in  the  nature  of  the  closure 
relationship.  They  justify  their  method  of  closure  by  relating  it 
to  an  expansion  in  the  coupling  parameters;  but  the  limitations  of 


their  technique  are  not  clear  and  a  method  of  improving  upon  their 

results  is  not  obvious. 
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The  procedure  followed  by  Vidal,  Cooper,  and  Smith  *   was  not 

as  inclusive  as  the  one  used  by  Capes  and  Voslamber,  but  it  does 
have  one  important  advantage:  the  expansion  in  the  density,  used 
in  order  to  get  an  expression  for  the  memory  operator,  is  well 
understood  and  can  be  related  to  a  diagrammatic  expansion  of  the 
self -energy  operator. 

In  this  dissertation  we  combine  the  advantages  of  the 
hierarchy  approach  of  Capes  and  Voslamber  with  those  of  the  weak 

coupling  limit.  In  the  remainder  of  this  chapter  we  will  develop 
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a  formally  exact  method  to  close  the  hierarchy.  *    In  Chapter  IV 

we  will  apply  the  weak  coupling  limit  to  the  formally  exact 

expression  for  D(a;t).  It  will  be  seen  in  Chapter  IV  that,  depending 

on  how  the  weak  coupling  limit  is  applied,  we  can  reproduce  several 

of  the  existing  theories  of  electron  broadening  as  approximations 

to  the  exact  theory.  Furthermore,  we  then  generate  our  own 

approximation  procedure  which  enables  us  to  develop  a  more  systematic 

theory. 

III-B    Formal  Closure 

First  Equation  of  the  Hierarchy 

Rather  than  immediately  employ  an  approximate  closure 
relationship,  we  will  proceed  formally  to  obtain  an  exact  relationship 
by  observing  that  C(a,l,,,#sjt)  represents  a  linear  map  of  an 
atomic  function  onto  a  space  containing  functions  of  electron 
coordinates  as  well  as  atomic  coordinates} 


•jLV 


tf(a,l***sjt)=U(a,l"'s;t)p(a)d, 

ff(a,l—ejw)=«(a,l—«ju>)p(a)J,  (III-6-1) 

U(a,l«««s;t)=Txs+1.#.N  Ve1     p^lp^Ca), 

0C^l—«i«K»t41.#-II{w"  ^  PaN}P"1(a).  (m-e-e) 

For  the  specific  case  where  the  trace  is  taken  over  all  N  electrons 
we  get 

I>(a$t)*U(ait)p(a)i 

5(a;co)=U(ajo))p(a)i  (lII-e-3) 

U(a|t)^x1.#,N{VNe"iLtp^}p"1(a), 

The  first  step  in  effecting  our  closure  will  be  to  eliminate  p(a)a 
in  III-B-1  in  favor  of  D(ajt)  in  equation  III-6-3.  Hence,  assuming 
that  an  inverse  exists  for  U(ajt)  and  G(ajco),  the  functionals 
which  result  are 

&(a,l—a|t|^(ast))=U(a-l—sjt)lT'1(a|t)^(a}t),  (III-B-6) 

or  the  Laplace  transformed  version, 

D(a,l'»«s;w|D(a;a)))=U(a,l«'«sjoj)U~1(ajw)D(a;to).  (III-B-6) 

U(a,l,,#sjt)  and  UCa^l^^sjoOCr^ajw)  will  in  general  be  extremely 
complicated  operators.     Equation  III-B-6  gives  us  a  formal  method 


of  closing  the  hierarchy  of  equations  at  any  level  (any  value  of  s). 

The  first  equation  of  the  hierarchy,  III-A-10,  may  therefore 
be  written 

(|Y+iLa):B(ajt)=-inTr1{Lea(l)U(a,ljt)U^(ajt)}S(ajt).    (III-B-7) 
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This  can  be  cast  in  a  different  form, 


(3£.+iLa)D(ait)=-i 


'  dVx(t'-t)D(ajt)  (III-B-8) 

0 


by  introducing  a  collision  operator,     x(t'-t),  which  is  non-local 
in  time.    This  is  a  particularly  easy  form  to  Laplace  transform. 
Inspection  of  equation  III-£-£  shows  that  as  t->  0,  6(a,l»«»sjt) 
approaches  a  time  independent  functional  of  S(ajt) 

£im£(a,l»-*sjt|D(a;t)=£(a,l*»'S;t=Op(a;t)) 
t+0 

=U(a,l«"Sjt=0)u"  (ajt=0)D(a;t).  (III-B-9) 

This  shows  thatx(t'-t)  has  a  singular  contribution  at  f^t. 
Extracting  this  part  explicitly  fromx(t'-t)  yields 

(§£  +iLa)D(a;t)=-iBD(ast)-iJ  dt'M(t--t)S(ajt')  (lII-B-10) 

where  B  is  time  independent.    The  operator  M(t'-t)  is  now  non- 
singular  and  the  integral  in  which  it  appears  vanishes  as  t>  0. 
If  equation  III-e-10  is  Laplace  transformed  we  find  that 


( w-L  )D(  a;  w)=iD(  a; t=0) +1~4(  u)6"(  a j  to) ^C a;  w) 

3 


(in-e-ii) 
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This  result  displays  the  sane  functional  form  obtained  by  VCS  where 
M(to)  together  with  3  plays  the  role  of  the  VCS  memory  operator. 
It  can  be  shown  that  B  vanishes  when  electron-atom  interactions 
are  neglected  from  the  density  matrix.     Consistent  with  the 
separation  of  the  singular  part  of  the  collision  operator,  x(t'-t), 
from  the  non-singular  part,  we  write  "5(a>l««,sjt|D(ajt))  as  the  sum 
of  its  short  time  limit,  equation  III-£-9,  and  a  time  dependent 
remainder  which  vanishes  at  t=0: 

D(a,l»»'Sjt|D(a;t))=U(a,l«««sjt=0)u""1(ajt=0)r?(a;t)+  ^a,l"«sjt)(lH-B-12) 
where 

?(a,l"#$it=0)=0  „  (III-B-13) 

It  can  be  seen  that  ?(a,ljt)  is  related  to  the  non-singular  part 
of  the  collision  operator,  x(t'-t).    Equation  III-B-9,  which  is 
sometimes  called  the  short  time  limit,  is  discussed  in  Section  C 
of  this  chapter. 

We  now  make  a  few  observations  about  tf(a,l»«»s;t=0|t>(a;t)) 
and  ?(a,l»»»sjt),  given  by  equations  III-B-9  and  III-B-12. 
From  the  definition  of  U(a,l,#,s;t)  given  by  III-fc-2  we  have 

U(a,l...s;t=0)^)(a,l*"s)p'"1(a)j  U(a,t=0)=l 
U(a,l-5jt=0)u"1(a,t=0)^(a,l«"s)p"1(a)  .  (III-B-14) 

Substituting  III-B-14  into  III-B-12  and  taking  the  Laplace  transform 
gives 

rj(a,l"9jw|D(a;w)=p(a,l»"s)p~  (a)D(ajw)+p(a,l"»sjw).        (lll-e-15) 
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The  closure  relationship,  equation  111^-6  in  conjunction  with 
equation  IIl-e-15  enables  us  to  write: 

?(a.l"-S5w)=K(a,l«"SSio)5(ajw)  (III-B-16) 

K(a,l'«-s;co)=U(a,l--»sja3)u"*1(a;a))-p(a,l«-s)p"1(a).       (lII-fi-17) 

The  functional  form  of  the  operators,  B  and  M(wl  can  now  be 
exhibited  explicitly.     Considering  the  Laplace  transform  of  the 
first  equation  of  the  hierarchy,  equation  III-A-10, 

(co-L  )3(ajw)=hTr  {L     (l)5(a,l;w)}=ip(a)3 
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and  by  using  equations  III-6-15  and  III-C-16  we  arrive  at 
(u-L  ^(ajujHiTr^L    (Dpta^Dp^CaJ^aja)) 

-nTrJL    (l)K(a,l;w)}fj(a;w)=iP (a)c[  .  (III-B-18) 

This  can  be  compared  to  equation  III-B-11  to  get 

0(3^)= i~T P(a)d~ 

w-L  -B-M(cj) 
a 

B^Tr.U    (DpUDp^U)} 
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M(a))=nTr(Lea(l)B(a,l;a))}  .  (lII-B-19) 

The  operator,  B,  as  it  is  defined  in  equation  III-e-19,  is  a  simple 
functon  of  well  defined  operators  and  can  be  calculated.  The  operator 
M(co),  on  the  other  hand,  contains  the  formal  operator,  K(a,ljuj), 
which  is  obtained  from  the  next  equation  of  the  hierarchy. 


Second  Equation  of  the  Hierarchy 

The  approach  used  to  cast  the  first  equation  of  the  hierarchy 
in  the  form  of  a  linear  kinetic  equation  is  well  known.  However, 
it  has  been  shown  that  if  we  close  the  second  equation  of  the 
hierarchy  in  the  same  manner  as  the  first,  a  useful  result, 
which  may  be  easily  and  systematically  approximated,  is  obtained,* 
The  approximation  technique  used  will  be  discussed  in  detail  in 
Chapter  IV. 

From  equation  III-A-10  we  have  the  second  equation  of  the 
hierarchy, 
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(9_+iL(l))D(a,l;t)=^inTr2{(LQo(2)+L^)D(a,l,2jt)}. 

at 


ea 


ee 


In  order  to  put  this  into  a  form  where  k(a,lsw)  appears  explicitly, 
we  first  use  equation  III-B-15: 


S(a,l;t)=p(a,l)p"  (a)D(a;t)+P(a,ljt) 


B(a,l,25t)=p(a,l,2)p~  (a)S(ajt)+P(a,l,2st). 


(III-B-20) 


If  these  relations  are  substituted  into  equation  III-B-20  the 
resulting  expression  is 

(1_  +iL(l))P(a,l)p"1(a)D(a;t)+  (£_  +iL(l))p(a,ljt) 
=-inTr2{(Lea(2)+L^e1)P(a,l,2)p**1(a)}D(a;t) 


-inTr9{(L     (2)+L^)?(a,l,2;t)}. 

fc  6  3  6S 


(III-B-21) 


This  equation  is  developed  further  in  Appendix  A  where  the 
properties  of  the  equilibrium  hierarchy  for  an  N+l  particle  system 
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are  used.    The  resulting  kinetic  equation  is 

G_+iL(l))P(a,l;t)=-iGL  a(l)p"1(a)5(ajt)+ip(a,l)x 

ST  ea 


xp  (a)nTr1{Lea(l)P(a,lst)}-inTr2{(Lea(2)+L^)t*(a,l,2st)}(III^^2) 
where  the  operator  GL  (1)  is  given  by 

6  3 

GLea(Dp"1(a)3(ajt)=[p(aJll)-np(a,l)p"1(a)nTr2{p(a,2)P21} 

Wr2{p(a.l,2)P2i}]Lea(l)p-1(a)5(ajt)  (IIM-23) 

and  where  P. .  is  the  permutation  operator,  defined  in  Appendix  A. 
Taking  the  transform  of  equation  IIl«*-22  we  have 

(u)-L(l))?(a,lja3)+p(a-l)p"1(a)nTr1{L  (l)P(a,ljo))} 
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^Tr2{(Lea(2)+Le2e1)?(a,l,2;a))}^Lea(l)5(a3a))  (XII-M4) 


where  we  have  used  equation  III-B-13.  If  we  use  the  definition  that 
P(a,l,»#s;u))=K(a,l,»»sja))6(ajw)  then  equation  III-B-24  becomes 

(w-L(l))fea,l;a>)-»p(a.l)p"1(a)nTr1{L  (l)K(a,lju)> 


^Tr2{(Lea(2)+Lee)^a'1'2^)=GLea(l)  •  (IH-B-B) 

Thus  we  have  a  kinetic  equation  for  K(a,lju))  in  terms  of  K(a,l«2;<o  ) 

and  a  source  term,  GL    (l).    We  now  make  the  formal  definition 

ea 

V(a,liw)K(a,l;a))=hTr2{(L    (2)+L21)K(aJ1.2;u)> 
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-p(a*l)p""1(a)nTr1tt    U)K(a,lju>)}  (HI-e-26) 
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which  is  discussed  in  Appendix  B.    With  the  use  of  equation  III-B-26 


we  solve  equation  III-B-25  formally  to  get  a  convenient  form  for 
K(a,l;J 


K(a,l;u>)= rryz GL 

(o-Lu;-V(a,l;u) 


3 


(in-e-27) 


The  memory  operator  then  becomes 


MM-inTr^L     (1)         (1J  j  GL    (1)}, 


(III-B-28) 


It  is  now  possible  to  analyze  V(a,l;u))  in  exactly  the  same  manner 
that  we  analyzed  the  collision  operator,  x(t"-t),  that  is,  we 
separate  V(a,ljuj)  into  frequency  dependent  and  frequency  independent 
parts : 


V(a,lH=V(a,lja>=«>)+V  (a,ljto)  • 

c 


(lII-S-29) 


The  infinite  frequency  term,  VCa,^^),  is  analyzed  in  Appendix  B 
with  the  result  that 


V(a,ljaj=*>)K(a,l;w)=tfrr0{(L  (2)+L21)x 

2   ea    ee 


-l 


xX(a,l,2;a3=-)%"  (a,l;a3=co)}K(a,l;w)-p(a,l)p'"1(a)Tr1{L  (l)K(a,l;(o) } 


K(a,l»«»Sjto)=%(a,l»«»sjio)L  (l). 

63 


(III-B-30) 


The  collisional  part  of  111^3-29  can  best  by  analyzed  by  continuing 
to  the  next  equation  of  the  hierarchy. 
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Third  Equation  of  the  Hierarchy 

We  already  have  expressions  for  M(cj)  and  J(Ue)  in  terms  of 

V(a,ljw=:=°)  which  are  formally  more  inclusive  than  any  previously 

derived.  In  order  to  fully  understand  what  part  of  V(a,lju>)  is 

included  in  its  short  time  limit,  V(a,lju}ISo),  we  must  understand 

what  is  excluded  from  it;  that  is,  we  must  look  at  V  (a,ljio).  To 

c 

do  this  we  will  require  the  third  equation  of  the  hierarchy, 

(~  +iL(2)J6(aJ,lJ2;t)-inTr3{  (Lea(3)+L**+l£* )B(a*l,2«3|t)}  .  (III-A-IO) 
From  Appendix  B,  equations  B-9  and  B-10,  we  know  that: 
k(a,l'"s;a))=£(a,l"»s;w)L  (l)p  1(a) 

©a 

k(a,l*"sju))=*;(a,l'"s;u))£  X (a,ljco)K(a,l;o>).  (III-B-31) 

Then  in  equation  B-ll  we  defined  the  operator  V(a,l;w): 

V(a,ljco)K(a,l;co)=nTr  {(L  a(2)+L2Mx(a,l,2;w)^(a,l;a))}K(a,ljw) 


-p(a,l)p   !(a)NTr  |Lea(2)P2i}k(a,l5a)).  (^-l?.; 

We  now  want  to  divide  V(a,l;co)  into  frequency  dependent  and 
independent  parts,  in  a  manner  suggestive  of  the  separation  of  the 
collision  operator  in  equation  III-B-llo  In  order  to  do  this  we 
consider  the  inverse  transform  of  K(a,l***sjcj)  and  write  it  as  a 
functional  of  the  inverse  transform  of  K(a,lju>), 

K(a,l*"sjt)=K(a,l"'Sjt|K(a,ljt)).  (III-B-32) 

Following  a  procedure  analogous  to  that  used  to  derive  equation 
III-B-12  we  now  separate  this  last  expression  into  a  short  time 
functional  and  a  time  dependent  remainder  term: 


K(a,l"»sjt|  K(a«ljt))eK(a»l**««|teO|  K(a,l;t))+X(a,l—8;t).  (III-B-33) 

The  Laplace  transform  of  equation  III-3-33  yields 

k(aJl»««s;a3|K(a,l;w))=K(a,l-««s;t=0|K(a,l;a)))+X(a,l-"sju).  (III-B-34) 

Hence,  comparing  this  result  with  equation  B-13,  we  observe  that 
the  infinite  frequency  limit  is  related  to  the  short  time  limit: 

~K(a,l»«*sju|K(a,ljaj))=?<:(a,l«"S;co=*»)X,,,i(a,l^=o°)K(a,l^) 

+X(a,l'«*s;to)  (III-B-35) 

X(a,l»«»s;co)=x(a,l,,*s5a3)K(a,l;a))=[X(a,l«"sju))>£:",1(a,l;a)) 

-^(a,l«««S50)=»)  Aft(a,l«*«sja)^o)]K(a,l}cJ).      (lII-S-36) 

The  similarity  of  the  above  with  equations  IIIH3-12  and  III«£-17 

suggests  that  the  approach  used  to  determine  the  functional  form 

of  V  (a,ljcu)  from  the  third  equation  of  the  hierarchy  will  follow 
c 

that  used  in  Section  III-*  to  find  the  functional  form  of  M(co) 
from  the  second  equation  of  the  hierarchy. 

With  the  exception  of  an  intermediate  step,  the  algebra  involved 
in  converting  the  third  equation  of  the  hierarchy  into  one  for 
x(a,l,2;to),  defined  by  equation  III-*-35,  is  similar  to  that 
used  to  convert  the  second  equation  of  the  hierarchy  from  an 
equation  for  l)(a,ljt)  to  one  for  K(a,l}a,).32*33  The  result  of 
this  conversion,  determined  in  Appendix  C,  is 

(w-L(2bx^a,l,2;to)^(a,l,2;uoo)%  Ha.liw8*.)* 

xInTr2^Lea(2)+I?1eS(a*1.2ja))]-nTr3{(Lea(3)+L^+ 

+Lee)X(a*1'2;aj)}=S(a*1'2)  (III-B-37) 

where  S(a,l,2)  is  a  complicated,  frequency  independent,  source  term. 
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If  we  further  define  an  operator,  W(a,l,2;m),  as  follows, 
W(a,l»2ju)x(a,l,2ju)=X(a,l*2j<o=»)^  ^a.lju^.Jx 

"nTr2{(Lea(2)+Lee)X(a'1'2j(ij)}^Tl3{(Lea(3)+Lee  + 

+Le|)x(a*l*2*3jw)  ,  (III-B-38) 

we  can  formally  solve  for    x(a,l,2;aj)  and  V  (a,ljto)  in  terms  of 

c 

S(a,l,2)  and  W(a,l,2ju)): 

x(a,l,2;co)= rrr S(a,l,2), 

a)-L^;-W(a,l,2jco) 

Vc(a,ljo3)=-inTr2{(J2)+L^)  j^y- S(a,l,2)},      (III-B-39) 

oi-Lv     -W(a,l,2;u) 

where  we  recognize  that  V  (a,ljw)  is  still  formally  exact  with  all 

c 

the  complicated  N-body  effects  hidden  in  the  operator  W(a,l,2;co). 


HI-C    Short  Time  Limit 

In  the  next  chapter  we  will  consider  some  possible  weak 
coupling  limits  to  the  memory  operator,  M(co).    Before  we  do  this, 
however,  it  will  be  useful  to  consider  an  alternative  approximation 
method.     In  equation  III-B-10  we  rewrote  the  first  equation  of  the 
hierarchy  in  the  form: 

rt 


(  hr+il  )$(ajt)=-iBD(a;t)-i 


dt'M(t'~t)D(a;t'),  (III-B-10) 

0 


where  B  represented  the  singular,  time  independent  part  of  the  collision 
operator,  x(t'-t),  and  M(t'-t)  represented  its  time  dependent, 
nonsingular  part.  From  the  Laplace  transform  of  equation 
III-B-10  we  got 
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ft->-„-l.-MUP,*'  (IM"10) 
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B  was  called  the  short  time  limit  of  xCt-*-^)  and  M(oo)  vanished  in 
the  infinite  frequency  limit.  Hence,  the  memory  operator  was  divided 
into  a  frequency  dependent  and  a  frequency  independent  part.  We 

will  see  in  this  section  that  it  is  useful  to  apply  the  analogous 
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separation  to  higher  order  equations  in  the  hierarchy. 

Equation  III-4J-11,  above,  is  exact  and  has  the  same  form  as 

the  expressions  for  M(w)  and  V  (a,l;u)  which  are  also  exact: 

c 

MfaJs-inTrjU    (1)  ,    i  ;      GL    (l)p   ^a)}       (III-B-28) 


c 

.21,  i 


S(a,l,2)}.  (III-B-39) 


Vr(a,l;w)=-inTr?{(L     (2)+L     )  rev-- ~ 

c  2      ea  ee     wLlZ;-W(a,l,2iw) 

Thus,  in  formally  closing  the  first  equation  of  the  hierarchy,  we 

have  put  all  the  effects  of  the  N  electrons  into  the  atomic 

operators,  B  and  M(w).  Similarly,  in  formally  closing  the  second 

equation  of  the  hierarchy,  we  have  cast  the  evaluation  of  M(w) 

into  the  form  of  an  effective  two  body  problem  (the  atom  and  one 

electron)  where  the  operators,  V(a,l;o)=°°)  and  V  (a,l;co)  contain 

c 

the  effects  of  N-l  electrons.     Thus,  in  both  the  expression  for 
S(aju )  and  that  for  M(w)  the  many  body  effects  have  been  divided 
into  a  short  time  (infinite  frequency)  limit  and  a  frequency  dependent 
part  that  vanishes  at  t=0  ("=").     The  former  terms  contain  mean 
field  effects  while  the  latter  contain  collisional  effects. 

Next  we  consider  the  physical  significance  of  B  and  M(w)  as 
they  appear  in  the  effective  atomic  resolvent  which  governs  the 
time  development  of  the  operator  p   d,  equation  III-£-ll.     If 
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we  were  to  neglect  M(u)  then  all  electron-atom  effects  would 

be  contained  in  B  and  the  time  development  of  p  d  would  be  governed 

3 

by  the  mean  field  electron-atom  interactions.    The  analysis  of  the 

separation  of  V(a,lju))  into  V(a,lsar*»)  and  V"  (a,ljuj)  in  equation 

c 

III«£-28  is  very  similar.  If  we  were  to  neglect  V  (a,ljw),  the 

c 

time  development  represented  by  the  resolvent  operator, 
equation  III-e-28,  would  be  governed  by  the  exact  interaction  of 
the  atom  with  one  electron,  1/  ',  plus  the  mean  field  effects 
of  the  remaining  N-l  electrons  included  in  V(a,l;w=*>°).  Another 
way  of  stating  this  is  that  the  atom  is  perturbed  by  a  single 
electron  moving  in  a  static  background  due  to  the  other  N-l  electrons. 

This  discussion  suggests  a  possible  systematic  approximation 
method  for  calculating  the  line  shape  function,  l(w).  Keeping 
only  B  in  the  expression  for  I  (to)  is  the  crudest  approximation  and 
leads  only  to  a  slight  shift  in  the  location  of  the  line  center. 
Making  the  approximation  at  the  next  level,  that  is,  neglecting 
only  V   (a,lju)  from  the  denominator  of  the  expression  for  M(w) 
includes  the  interaction  of  the  atom  with  one  electron  exactly, 
together  with  the  mean  field  effects  of  the  other  N-l  electrons. 
It  will  be  observed  in  the  next  chapter  that  most  existing 
theories  of  electron  broadening  can  be  obtained  by  taking  some 
sort  of  weak  coupling  limit  of  the  result  of  this  approximation 
to  M(o>).  It  is  possible  to  go  further  and  keep  only  the  short 
time  limit  of  W(a,l,2ju>)  which  appears  in  the  denominator  of  the 
expression  for  V  (a,l;w),  but  since  retaining  only  V(a,l;uF°°) 
yields  an  expression  for  the  memory  operator  which  goes  beyond 


<IU 


most  existing  theories,  the  nature  of  V  (a«ljco)  will  not  be  emphasized 

c 

here. 

The  approximation  procedure  outlined  in  this  section  should 
be  regarded  as  a  possible  alternative  to  the  weak  coupling  limit. 
It  is  hoped  that  this  method  will  serve  as  a  starting  point  for  a 
fully  quantum  mechanical  unified  theory  applicable  to  a  high 
density  plasma. 


CHAPTER  IV 
WEAK  COUPLING  LIMITS 

IV-A  Introduction 


In  Section  III-C  we  discussed  an  exceedingly  useful 
procedure  for  obtaining  approximate  line  shape  functions  by  taking 
the  short  time  limits  of  the  different  collision  operators.  In 
this  chapter,  however,  we  will  consider  a  different  approximation 
procedure  which,  we  shall-  show,  parallels  more  closely  the 
results  of  existing  theories  of  line  broadening.  This  technique, 
called  the  weak  coupling  limit,  involves  a  perturbation  expansion 
in  some  sort  of  coupling  parameters;  the  impact  approximation 
mentioned  in  Chapter  II  is  related  to  this  technique.  In  this 
chapter  we  identify  possible  expansion  parameters  and  explore 
soma  expansion  techniques.  We  will  find  that  variation  in 
the  weak  coupling  methods  lead  to  different  expressions  for  the 
memory  operator,  some  of  which  we  will  relate  to  existing  theories, 
One  expansion  will  be  developed  that  leads  to  an  expression  for 
the  line  shape  that  contains  several  of  the  previously  developed 
theories  as  approximations,  and  which  can  be  systematically 
carried  further  than  any  of  them. 

In  Chapter  III  the  results  were  exact  within  the  limitations 
of  the  static  ion  approximation;  neither  the  classical  path  nor 
the  no-quenching  approximations  were  made.  However,  the  bulk  of 
the  literature  concerning  line  broadening  has  been  within  the 
framework  of  these  approximations,  and  they  will  also  be  assumed 
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in  the  remainder  of  this  dissertation.     Ch  the  other  hand  it 
will  be  shown  that  several  of  the  important  approximations  which 
have  been  introduced  in  an  ad  hoc  manner  by  many  authors,16*17*22*23 
the  impact  approximation,  the  neglect  of  electron«electron 
correlations,  and  the  neglect  of  electron-atom  correlations  from 
the  density  matrix,  will  follow  from  simple  expansions  of  the 
memory  operator  in  the  various  coupling  parameters. 

However,  a  difficulty  that  arises  in  the  coupling  constant 
expansion  as  applied  to  the  line  broadening  problem  in  the  existence 
of  two  intrinsically  different  types  of  interactions:     the  electron- 
electron  interaction  and  the  electron-atom  interaction.    The 

electron-electron  interaction  is  treated  extensively  in  the  plasma 

37  38 
theory  literature.     *        The  electron-atom  interaction,  however, 

will  contain  atomic  operators  which  are  unrelated  to  the  fundamental 

lengths  of  the  usual  plasma  problem.    This  suggests  that  the 

coupling  constant  we  use  for  the  electron-atom  interaction 

should  be  independent  of  the  one  we  use  for  the  electron-electron 

interaction.     In  plasma  line  broadening,  the  concept  of  truncating 

a  hierarchy;  of  kinetic  equations  by  expanding  simultaneously  in 

two  independent  parameters  was  carried  out  first  by  Capes  and 

27 

Voslamber.        While  also  using  a  two  parameter  expansion,  we 

will  employ  a  different  electron-electron  coupling  parameter 

and  will  apply  the  weak  coupling  limit  in  a  sojoswhat  different  way. 

We  now  discuss  possible  expansions  and  expansion  parameters 
further.  The  most  straightforward  method  is  to  expand  in  powers 
Of  the  coupling  constants,   X       and  A     ,  defined  by  the  relations, 

66  63 
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Vee(^eeVee(j)'     Wj)=Vee(^  (lV-A~2) 

This  is  the  approach  used  by  Capes  and  Voslamber  and  we  will  later 
show  that  our  technique  for  approximating  the  kinetic  equations 
gives  a  result  very  much  like  theirs  if  we  go  to  the  same  order 
in  the  two  coupling  constants  as  they  did.     The  expansion  in  the 
electron-atom  coupling  constant  is  a  good  one,  as  seen  in 
Appendix  D,  and  we  will  use  it  in  this  chapter.     On  the  other  hand 
it  has  been  shown  in  the  literature  that  for  large  ranges  of 
plasma  temperature  and  density  where  the  electron-electron  coupling 

constant  may  not  be  a  valid  expansion  parameter,  an  expansion  in 

37 

the  plasma  parameter  is  the  better -choice. 

The  procedure  involved  in  expanding  in  the  coupling 
constant  for  an  interaction  is  a  simple  one  which  orders  contributions 
to  quantities  being  expanded  in  powers  of  that  constant.  The 
concept  involved  in  the  plasma  parameter  expansion,  however,  is 
slightly  more  subtlej  before  the  coupling  parameter  is  ever 
identified,  all  expressions  are  first  scaled  to  an  appropriate 
set  of  units.  The  natural  units  with  which  to  scale  the  plasma 
problem  are 

!=(47rne2/a)§'i     LpHw/^neZ)? 
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x*x'=  f-  j  or*.)'  =  ~  ;     v*v=  r-*-  .  (IV-A-3) 

T)  wp  Vp 

As  shown  in  Appendix  D,  where  the  equations  involved  in  the  line 
broadening  problem  are  scaled  to  these  lengths,  the  appropriate 


expansion  parameter  is  the  plasma  parameter: 

where  N  represents  the  number  of  particles  in  the  Debye  sphere. 

It  is  also  observed  that  if  the  kinetic  equations  of  line  broadening 

are  scaled,  the  electron-atom  coupling  parameter  becomes 


From  Chapter  III  the  line  shape  due  to  electron  broadening 


Xea=  — :  (IV-A-5) 

V 


is 


^E>=KTra{3  ^nbra-p(a)J}  (IV"A-6) 

a 
where  the  factors,  p  ,  B,  and  M(a>),  will  in  general  contain  terms 

3 

which  include  all  orders  in  all  of  the  coupling  parameters. 
Expanding  p(a),  which  appears  in  the  numerator  to  lowest  order 
in  X   and  x   is  equivalent  to  replacing  it  by  e"°^  /Tr{e   }, 
where  H  is  the  unperturbed  atomic  Hamiltonian.  The  resolvent 
operator    »  -b-\T(  )  *s  more  difficult  to  approximate,  however. 
A  straightforward  perturbation  expansion  which  is  first  order  in 
any  of  the  coupling  constants,  yields 
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=  -V  jl-KB+MU))11'  -V  ]  (IV-A-7) 
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where  (B+M(to))  represents  the  lowest  order,  nonvanishing,  term 

in  a  perturbation  expansion  of  B+M(u).     Near  the  line  center,  where 
M((jj)  approaches  M(0),  the  effective  atomic  resolvent  appearing  in 
the  line  shape  function  becomes     — rTSTTToT  »     Since       *r 

(0-«L  -4$— n\.v)  bi—L 
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is  not  necessarily  small,  an  expansion  of  the  form  shown  in 
equation  IV-A-7  will  not  be  a  good  one.     *        Thus  B  and  M(u), 
which  may  be  viewed  as  "width  and  shift"  operators  for  the  line 
shape,  should  be  retained  in  the  denominator;  it  is  then  possible 
to  expand  B  and  M(oj)  themselves  in  the  various  coupling  parameters. 
In  the  remainder  of  this  chapter  we  will  discuss  the  possible 
methods  of  expanding  them. 

From  Chapter  III  we  recall  that  the  operators  B  and  M(u) 
are  given  by 

B=hTr  {L .     (l)p(a,l)p    l(a)  (lII-B-25) 

1      6  3 

M(J=-inTr  {L .     (1)  jrr^z GL     (l)p~1(a).  (III-B-28) 

1     ea  T vx/  ,,/     ,      \       ea      H 

The  expansion  of  B  in  the  various  coupling  parameters  is  a  straight- 
forward problem,  but  M(u),  which  contains  the  formally  exact 
operator,  V(a,lju)),  is  difficult  to  deal  with.     It  is  the 
denominator  appearing  in  equation  III-B-28  that  causes  the 
difficulty.     Again,  we  could  expand  the  resolvent  in  the  manner 
suggested  by  equation  IV-A-7,  but  this  expansion  would  also  be 
invalid  in  the  line  center.     Therefore,  we  examine  the  operators 
in  the  denominator,  lr     +V(a,lju),  in  the  various  weak  coupling 
limits.    The  results  of  these  approximations  to  the  denominator  of 
the  expression  for  MCo )  are  displayed  in  Table  1  and  are  discussed 
in  detail  by  the  remaining  sections  of  this  chapter. 
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TABLE  1.  Some  Weak  Coupling  Limits  for  the  Denominator  of  M^) 


Some  Weak  Coupling  Limits.. 
for  the  Denominator  of  M(oj) 


RESULT 


*ea%e=° 


second  order  theory 
no  shielding 


X  =0;  A  =0 
ea 


second  order  theory 

"random  phase  approximation" 


Xea=lj     Xee=0 


unified  theory  of  VCS 


K a=1S  AaS0  for  G; 

X    =1  for  denominator  of  M(u)) 
ee 


unified  theory  of  Capes 
and  Voslamber 


Xea=l;  A  =0 


unified  theory 

"random  phase  approximation" 


W-B     Second  Order  Tjjggrjgg 
The  simplest  class  of  approximations  for  G  and  M^)  in 
equation  IIIH3-28  are  those  which  retain  only  the  lowest  order 
nonvanishing  terms  in  an  electron-atom  coupling  constant 
expansion.     Second  order  theories  are  so  named  because  by  taking 
L      "*V(a*l;u))  and  G  to  zeroth  order  in  x       we  restrict  the  memory 

6  3 

operator,  MGd),  to  second  order  in  the  electron-atom  coupling 

parameter.     Even  here,  however,  there  is  some  latitude  in 

approximating  the  electron-electron  interactions;  in  this  section 

we  will  consider  the  cases  where  Xrt_=n  and  A=n» 

ee 

As  the  first  approximation  to  both  G  and  1/     +V(a,l;u)),  we 

consider  X     =0  and  X     =0.     From  equation  III-E-8  we  have 
ee  ea 

I/1  ~L  -iv  *y   +L     (l)  which  to  zeroth  order  in  both  coupling  constants 
a      i   vi     ea 

becomes  V   '-iv  *V  .     The  operator,  GL _   (l)p'"1(a),  can  be  simplified 
li  ea 

if  we  realize  that  to  lowest  order  in  \     ,  p (a,l)=p (a)p (l)  and  to 

6  o 

lowest  order  in  X     ,  p  (l,2)=p  (l)p(2).     Hence  from  equation  III-G-23 

Uv 

GL    (l)p-1(a)=fn  (vJp(a)L .    (Dp-^a).  (IV-B-1) 

Since  all  contributions  to  V(a,l;(o)  contain  at  least  one  factor  of 
L1J  or  LI  (j),  thus  this  operator  will  not  contribute  in  this 
approximation  and  R(o))  becomes: 

M(w)s-infd(l)L  (1)  * f  (v  )p(a)L  (Dp-^a).    (IV-B-2) 

J     ea    a)-La+i^1^i   °   1      ea 

The  quantum  mechanical  analog  of  this  is  essentially  the  expression 

18 
that  was  calculated  by  Smith  and  Hooper. 

Another  result  is  obtained  by  taking  G  and  lr     +V(a,lju) 

to  zeroth  order  in  the  plasma  parameter  rather  than  to  zeroth  order 
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41-43  (l)  -*■    * 

in  X     •  "  Again,  in  this  approximation  we  have  Lx   '=L  -iv  #v  j 

but  G  and  V(a,l5cj)  will  both  be  more  complicated  than  they  were 

for  the  case  where   X    =0.     In  equation  D-12  of  Appendix  D  the 

operator  GL    (l)p~  (a)  is  scaled  to  dimensionless  coordinates; 

G  "t"  P     U)=  ^V~  [p(a)9(a'lheaLea(rr)+  A  fc'PW" 


P 


*(g(a,l,2)-g(a,l)g(a,l))Lea(r2')  Jp^a).  (IV-6-3) 


If  we  now  observe  that,  for  X    =0,  g(a*l,,,s)=g(l»««s)  and  if  we 
recall  that  H(r  r  ),  defined  by  G(r  r  )=l+h(r  r  )  is  the  pair 
distribution  function  which,  when  scaled  to  dimensionless  coordinates, 
is  proportional  to  A,  then  equation  D-12  may  be  written  in  the  form, 

G  i*iii  p-l{a)=  *J*£  trf.)^ttfl-)*|tfa-p(.)htfI*I1^IWitfl1  ]: 

p  p   "  ... 

«p",(a)=^"p(a))bal«a(?1-)p-1(.) 

&<VKa<*rW*>(l>V>W*a'>-  (IV"B-4) 

L    (r  ')  is  an  electron-atom  interaction  which  is  statically  shielded 
ea     l 

by  the  electron-electron  pair  correlation  function.     In  Appendix  E 
it  is  shown  that  K(a,l,2;ar^°|K(a,ljcd))  to  zeroth  order  in  A      is 
given  by 

K(a,l,2jw=^|k(a,ljw))=f  (v  )f  (v  )[(1+P     )(l+h(l2)) 

1  0       10      2  21 


Hence  from  equation  E-16 


d(3)f o (v2 )g(123)P3   ] (l~C)f~i (vj )R(a,ljJ. 


V(a,l;to=^)=-inv  "v  fA(v)6f-1(v  ) 

110       1  o  1 
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where 


C=N 


d(2)fQ(v2)D(l2)P2i; 
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C(12)  is  the  direct  correlation  function.  When  this,  along  with 
equation  IV-B-3,  is  substituted  into  the  expression  for  the  memory 
operator  we  find  that 


MLJLsZLka. 


^r^V*! '>*<»--  t  ^VV^W*  o(vr)Cf-  W) 

xp(a)fft(v/)L*  (v/)p"1(a).         (XV-B-5) 

However,  the  form  of  equation  IV-B-5  is  not  simple  to  evaluate  since 

37 
a  modified  Vlasov  operator,      equation  E-16,  appears  in  the  denominator. 

We  have  shown  in  Appendix  F  that  any  expression  having  the  general 

form  of  equation  IV-B-6  can  be  cast  in  the  form  of  a  dynamically 

shielded  electron-atom  interaction;  thus, 


»i 


"*K 


£LJl=  -isa 


u_ 


dv  dv 

l     i     ea*   l 


L     (r     Alw  -  *"**•        i     1  0     1 


u 


.  ifciff1-)p-,W. 


(IV-B-6) 


where,  in  this  case,  the  dynamically  shielded  interaction  is  given 
exactly  by  equation  F-25: 


Lea      I27F 


dt 


(0)/*      T 


-«ik»r 


) 


where     (kjw«L  )  is  a  frequency  dependent  dielectric  function 


e(0)(W>[X+ 

3 


d+    __i ?'fc(k)fn(v)].         (IV-B-7) 


a     i     l 

If  the  static  shielding  defined  by  equation  IV«*-4  is  substituted 
for  the  dynamic  shielding  of  equation  IV-e-6,  we  finally  arrive 


DU 


at  a  second  order  theory,  including  electron  correlations,  similar 
to  that  studied  by  several  authors.19*20*41"43 


IV-C  Unified,  Theories  -  Expansion  in  \ 

&?ret-h  Order.  ?JUee 

Another  useful  result  is  obtained  by  talcing  the  denominator 

of  the  effective  two  particle  resolvent  operator  in  M^)  to  first 

order  in  A   and  to  zeroth  order  in  a  .  Again,  note  that  it  is 
ea  ee 

only  the  denominator  of  M(w)  that  we  are  expanding  to  first  order 
in  A   and  not  the  entire  function.  With  this  expansion  we 

S3 

arrive  at  an  expression  for  M(w)  v/hich  includes  all  orders  in  the 
interaction  of  the  atom  with  only  one  electron.  Such  an  expression 
for  M(oj)  yields  a  line  shape  function  that  is  valid  in  the  line 
wings  as  well  as  in  the  line  center;  hence,  it  is  called  a  unified 
theory.  It  should  be  observed  that  the  result  obtained  in  this 
approximation  will  entirely  neglect  electron  correlations. 

In  equation  IV-B-1  it  was  determined  that,  zeroth  order  in 
Aeg  and  Aee,  G  is  given  by 

GL  (l)p-1(a)=f  (v)p(a)L  JUp^U).  (iV-C-l) 

ea         u  1     ca 

In  this  approximation  the  operator  1/  +V(a,ljaJ)  is  just  as  simple 
to  evaluate.  L  is  retained  in  full,  while  V(a,lja))=;V(a,l;a)::^)  + 
+V  (a,l;w)»  given  by  equations  B-13  and  111-43-39,  will  vanish  to 

G 

zeroth  order  in  \       and  \     .  The  resulting  form  for  the  memory 
operator  is 

MU)=-infd(l)Lea(l)  ^  ^^  (l)   fo(vi)p(a)Lea(l)p-Ha).   (IV-02) 

J  3   6      6  3 


The  above  will  be  recognized  as  the  same  result  that  was  obtained 

22  23 
by  VCS,  whose  theory  was  outlined  in  Section  II-B.  *    It  will 

be  recalled  that  their  theory  used  a  density  expansion  of  the 

expression  for  the  memory  operator.  The  advantage  of  their  result 

was  that  it  included  all  orders  in  an  expansion  of  the  interaction 

of  the  atom  with  only  one  electron;  while  its  primary  disadvantage 

was  the  neglect  of  the  correlation  effects  of  the  other  N-l 

electrons.  As  mentioned  in  Section  II-B*  VCS  partially  compensate 

for  this  neglect  by  including  curoffs  in  the  spatial  integrations. 

In  the  remainder  of  this  chapter  we  will  show  that  the  mean  field 

effects  due  to  the  N-l  additional  electrons  can  be  included 

systematically  via  an  expansion  in  either  X      or  A. 

First  .teter  in  *ee 

We  have  just  seen  the  expression  that  is  obtained  for  the 

memory  operator  by  taking  its  denominator  to  first  order  in 

X   and  zeroth  order  in  X  .  We  next  approximate  1/  ^(a^u)) 
ea  ee 

to  first  order  in  X   or  in  X  . 

ea      ee 

To  first  order  in  X  ,  1/  '  is  retained  in  full;  but  V(a,lj") 

ea 

is  complicated  and  hence  must  be  approximated.    The  mean  field 
part  of  this  operator,  given  explictly  in  Appendix  D,  equation  D-15; 

rd(2)(Lea(2)+Lee){p(a)f0(v1)f0(v2)[g(a,l,2)(l+P2l),>= 


V(a,ljuf=°°)=n 


d(3)f  (v  )g(a,l,2,3)P    ~g(a,l,2)n  d(2)f  n(v  Jg(a,2)P     ]> 

0      3  31  J  u      z  cl 

x{p(a)f  (v  )  [g(a,l)-*n  d(2)(g(a,l,2)-g(a,l)g{a,2))P2i]r1 


«-np(a)f  (v  )g(a,l) 


d(2)L (2)P     . 
ea         2i 


t>2 


We  have  purposely  not  written  this  equation  in  terms  of  the 
dimension less  coordinates  related  to  the  plasma  parameter  because 
we  want  to  expand  it  in  X   not  A.  With  the  requirement  that 

(its 

equation  D-15  be  taken  to  first  order  in  either  *   or  X    and 

ea    ea 

since  there  is  an  explicit  factor  of  either  L   or  L   appearing 
under  the  integral,  the  remainder  of  the  integrand  must  be 
zeroth  order  in  both  of  the  expansion  parameters.  Under  this 
restriction  g(a*l*#,s)=l,  and  equation  D-15  reduces  to 

V(a,l|u>^)=Jd(2)l2lf  (v  )P  =V(1)  (IV-C-3) 

J       W  0   1   21 

where  V(l)  is  the  usual  Vlasov  operator  and  P. .  is  the  permutation 

operator  defined  in  Appendix  B.     The  collisional  term,  V  (a,lju,), 

c 

is  given  by 

V  (a,l*.)=-in(d(2)(L  a(2)+L")  (9\  -  i S(a,l,2)l  (III-B-39) 

c        j     ea    ee  u_i>2'-w(a,l,2;u>) 

It  can  be  shown  that  S(a,l,2),  appearing  in  the  integrands  is  first 

order  in  either  >   or  i  I  thus  we  neglect  V  (a^lj^).  We  now 

have  a  result  for  the  denominator  of  M(co)  which  is  exactly  the 

same  as  that  derived  by  Capes  and  Voslamber.  If  we  further 

take  the  operator  GL  (l)p-1(a)  to  zeroth  order  in  a  (not  first 

order  in  A  )  the  resulting  operator  is  the  same  as  that  obtained 

in  Section  IV-*: 

L  (1)   f  (vr)  e 
G-^—  =-Vt  L'(%)p -1(a)  (IV-6-3) 

'°p     >T)    ea  1 


Sfil^&i***®!*!'^®*')-  (IV^"4) 
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Substitution  of  equation  IV-£-3,  together  with  equation  IV-fi-4  into 
the  general  expression  for  the  memory  operator*  equation  III-B-28, 
gives, 

MU)=-injd(l)Lea(l)    (il  ■  ■  f0(v1)P(a)L;ap-1(a)<      {rl*>4) 

which  is  identical  to  the  result  obtained  by  Capes  and  Voslamber. 
If  we  now  apply  a  technique  very  similar  to  the  one  they  used,  and 
which  is  developed  in  Appendix  F,  this  equation  becomes 


M(u))=-in 


d(1)]£    Jw f  (v  )p(a)L^p-1(a)         (IV-C-5) 

a       li     ea 


where  the  generalized  dynamic  shielding  is  given  by 

Lja(J)«fdfat*Let(-«)D-1(^«'j„)e-iJ'». 

It  is  interesting  to  note  that  to  lowest  order  in  A   the  above 

ea 


\—\ 


expression  reduces  to 


I?  (?  )=[d£ **- 


a> 

«(6«^.)^  %r  fdv.v,  -^-^ - ,       (iv-c-6) 

1       a     l 
where  this  last  quantity  has  the  form  of  the  usual  plasma  dielectric 
function,  with  o>,  replaced  by  u,-L   . 

It  should  be  noted  here  that  the  mean  field  operator  appearing 
in  the  denominator  of  equation  IV-O-4,  and  therefore  the  dynamic 
shielding  of  equation  IV-C-6,  differs  slightly  from  that  appearing 
in  the  denominator  of  equation  IV-B-5.  In  that  second  order  theory 
we  took  the  denominator  to  zeroth  order  in  the  plasma  parameter 


~JH 


and  derived  the  generalized  Vlasov  operator,  equation  E-16;  in  this 
case  we  have  expanded  in  the  electron««lectron  coupling  strength. 

IV -0     lfrj.fp.ed  Jhe or v  -  Random,  Phase  Approximation 

In  Section  IV -A  we  suggested  that  for  large  regions  of 
temperature  and  density  an  expansion  in  the  plasma  parameter 
would  be  preferable  to  an  expansion  in  the  electron-electron 
coupling  strength.    We  will  see  in  this  section  that  an  expansion 
in  this  parameter,  coupled  with  an  expansion  in  x    *  is  a  distinct 

6  3 

improvement  over  theories  presented  in  previous- sections 5  the 
result  is  an  all  order,  or  unified,  theory  in  which  the  effects  of 
electron  correlations  are  systematically  accounted  for. 

Since  all  operators  in  this  section  will  be  expanded  in  the 
plasma  parameter  we  will  express  them  in  terms  of  the  dimensionless 
coordinates xde fined  in  Section  IV-^A.     From  Appendix  D,  equations  D-6, 
D«6,  and  D-7,  we  recall  that 

01  6"(aja))=i(a)'-  -*■-»—-  **™'"  )      p  d, 
P  w_       ml        ai  a 

P        P  P 


&-=nfd(l)p(a)f  (v  )g(a,l)  -" p^Ca), 


wp 


J  01  (^ 


"p 


L     (1)  .  L     (1) 


0) 

wp      wp 


Since  the  operator  B  will  involve  only  a  static  shift,  it  will  not 
be  emphasized  in  this  section;  instead  we  will  emphasize  a  study 

"p  L  (1) 

From  equation  D-12  we  can  express  G  ■*■   p  (a)  in 

"p 


» 


dimension less  units: 


(«)s^^p^)Xeag(a.l)Lea(;/)+^jd?2'p(a)g(a.l.2)Lea(?20 


-g(a,l)  f 

A 


df  -g(a,l)L .    (?  -), 

1  ea     i 


where  the  spatial  parts  of  the  reduced  distribution  functions  have 

GL  (1)  ml 
not  yet  been  scaled.  Since  we  are  requiring  that   aa    p  (a) 


WP 


be  taken  to  lowest  nonvanishing  order  in  the  electron-atom  coupling 
strength,  which  means  that  g(a«l*"a)«g(l***a)«  equation  D-12  becomes 
L    (1)  f  (v   -) 


p  (IV-D-1) 


P  P  ^  (IV-D-1) 

44 
It  has  been  shown   that  reduced  distribution  functions,  when 

expressed  in  dimensionless  units,  may  be  easily  expanded  in  A: 
g(r  ,r  )=l+h(r  ,r  )=1+Ah(r  ;r  ') 

12  12  12 

g(?i,r'2,?3)=l+h(?i?2)+h(??3)+h(?/3)+h(^r^3) 

=L+Ah(?  '£  0+Ah(?  -$  -)+Ah(r*  <$  -)+A2h(r  ^0i?  '). 

12  12  23  123 

(IV-D-2) 

This  result  may  be  substituted  into  equation  IV-D-1  with  the  result 
G  ^  ~'U)=  $£-  ^a)x.a[Lea^r>+K2^.'?2'>l„(?2')1'>"1<«) 

P  u 

This  is  the  same  as  the  result  obtained  in  Section  IV«£. 

We  must  now  evaluate  V(a,l;w)=^(a,l;io=°°)+V  (a,ljto)  to  first 

c 


yo 


order  in  the  electron-atom  coupling  strength  and  zeroth 

order  in  the  plasma  parameter.     It  will  be  convenient  to  separate 

this  approximation  to  V(a,ljto)  into  two  parts: 

V(a,lju))=V **  °(a,l; u>)+V  °*  °(a,lj u) .  (IV-D-4) 

VMle  V       (a,l;w)  is  zeroth  order  in  both   X      and   An  V  *   (a,l;w) 


is  first  order  in  A       and  zeroth  order  in  A.     In  Appendix  D 

we"  have  expressed  the  frequency  independent  contribution  to  V(a,l;) 

in  dimensionless  units: 

^^^V^W'.^^    -H*"  ^{p(a)fo(v/)fo(v/) 

P  J  ^2         X12  *i 

x[g(a,l,2)(l+P21)+^    d(3')f0(v30(g(a.l,2,3)-g(a,lJ,2)g(a,3))P31)]}« 
*{p(«)f-0(vl1  [g(a,l)+  ^  fd(2')f  0(v20(g(a,l,2)-g(a,l)g(a.2))P21]r1 

-%3(a.l)(dxYd?2\a(rY)P21 

■ 

where  L        =i(-""i  ~T")#  ""** *~  .    We  now  identify  the  terms  in  the 
ee  9af  r  12      3^ 

above  equation  that  are  zeroth  order  in  both  parameters.     The 
terms  having  integrals  containing  L     (j)  certainly  can't  be  zeroth 
order,  so  we  will  consider  only  those  terms  which  do  not 
explicitly  contain  this  factor; 


V  *  (a>ljw=°°)=i 


9Tj       r- 


■  •         9V! 


[g(l2)(l+P2x)+  \  {d(3')fo(v3')(g(l23)-l)P3i]}{f0(v/) 


[1- 


k 

A 


d(2')ffl(v/)(g(l2)-l)P21]} 


-1 


(IV-D-5) 


In  the  last  equation  we  have  used  the  fact  that  when  \    -*0, 
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g(a,l»««s)->g(l«»«s).    Equation  IV-D-5  will  be  recognized  as  the 
modified  Vlasov  operator,  which  in  Appendix  E,  equation  E-17,  is 
shown  to  reduce  to 

V#°*°(«,li«**')«-lnv1^1fo(vi)Cf^1(vi). 

We  must  now  look  at  V  (a,ljw),  given  in  dimensionless  coordinates 

c 

by  Appendix  D  (equation  D-17), 

V  (a,lju))  f  i  *  i  a 

"p  J  A      ea     2  ^      ri2         ^  , 

up       T> 
in  order  to  determine  whether  it  contributes,  to  zeroth  order  in 

both  parameters.  The  operator,  S(a,l,2),  given  in  Appendix  C 

by  equation  C-ll  can  be  shown  to  be  first  order  in  either  A   or  A. 

Since  the  free  streaming  part  of  the  effective  resolvent  in  equation 

D-17  is  zeroth  order  in  both  parameters,  we  see  that  there  will  be 

a  contribution  to  V  (a,ljco)  which  is  first  order  in  X   and 

C  0  3 

zeroth  order  in  A,  but  no  contribution  which  is  zeroth  order  in 
both  parameters.  Thus  we  have  identified  the  contribution  to 
V(a,lju>)  which  is  zeroth  order  in  both  parameters: 

V^'^lsJ^inv^yv^Cf;1^).  (W-D-6) 

The  contributions  to  V  *°(a,l;oj)  may  be  determined  by  examining 
equations  D-16  and  D-17  more  closely. 

The  operator  V1*0(a,ljt0)  is  still  very  complicated,  but 
there  is  a  class  of  terms,  each  having  a  very  simple  functional 
form,  that  will  be  seen  to  have  a  straightforward  physical 
interpretation.     Let  us  now  write  out  the  contributions  to  V(a,l;^) 
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which  contain  integrals  with  a  factor  of  L  (l)  in  the  integrand: 
V  (ajljur^)  t  X 

; — ■  *r<V  -f  Lea(?r}{fo(vr)f0(v/)[^i2)(i+p2i) 


+  J-  fd(3')f0(v3')(g(123M)P31]}{f0(vi') 


d(2')f0(r2')(g(l2)-l)P2i]} 


-1 


-  "**■  |d(2')L    (r   OP 


ea 


2         21 


(IV-D-7) 


where  we  have  demanded  that  the  above  be  no  more  than  first  order 

in  X     •    We  now  substitute  equation  IV-£>-2  into  equation  IV-D-7 
ea 

and  keep  terms  no  higher  than  zeroth  order  in  A: 


V_(a,l$u)=<x>) 


=X 


ea 


dr   -dv   'L    (r   -){f  (v   ')f  (v   ')  [(*•  +h(.     .  . 

2        2     ea      2  0      2        0      1  A  1      2 


rVrV))* 


x(l+P      )+ 
21 


dWi^^Khtf^^Mf^^MZ^^+MZf/J) 


xp  ]}(i^)rl(v  )-^ 

31J  0         1  A 


dr    'dv   'L     (r    ')P 

2       2    ea     2       21 


(IV-D-8) 


where  the  results  of  Appendix  E  have  been  utilized.     Making  use 
of  the  symmetry  properties  of  L     (j),  together  with  the  definition 
of  C,  the  above  equation  simplifies s 


V  (a,ljar°o) 
a 

*P 


=  X 


ea 

r 


^■^W^'h11"^  i)c"f"1(v  i» 


-i/ 


%af0(vi  >  ^^^(r^'Jh^^OCl-f  0(v2)Cf  0  (v2)P21 


+  X    f  (v    ')  dr    'dr*  1     (x    ')f  (v    ') 
ea    o     l  2       2    ea     2       o     2 


d?,u,%i.(»V/p(w  „<v  pa-oH*  p). 


(IV-D-9) 

From  equation  IV-O-8  we  observe  that  all  but  one  of  the  terms  contain 
a  permutation  operator  either  explicitly*  or  implicitly  through  the 
operator  C.    We  will  separate  out  these  terms  and  rewrite 
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equation  IV-D-9  in  the  abbreviated  form, 

V  **  °(a,ljoF«>)  t 

-9 =X      dx  "dv  a     (r')h(r   "r  Of  (v   ')+X(P,.)  ,_.  -  ... 

w  ea       2       2     ea     21202  lj  UV-4J-10; 

P         J 
where  X(P. .)  represents  the  contributions  that  contain  P.  ..  We  now 
combine  equation  IV-*-10  with  the  expressions  for  V  *  (a,l;u)}, 
V**°(a,l|w=*°)*  and  V**°(a,l;w): 

6  C 

Wp  J 

A(co)^*0(a,i.(^»)+v^°(a,ljw)+X(Pij).  (IV-O-II) 

In  equation  IV-C-11  we  have  explicitly  separated  out  the  simple 

multiplicative  frequency  independent  term,  and  also  the  term  which 

is  zeroth  order  in  both  expansion  parameters,  from    ■  AiUAiittf.  # 

Hence,  A(to;,  which  is  first  order  in   X    ,  contains  integral  operators, 

ea 

some  of  which  will  be  frequency  dependent.     We  now  explicitly  write 
out  the  memory  operator,  with  V(a,l;w)  given  by  equation  IV-O-11: 


a*-*  &  fd(i')LMc*/)i(^-^  *v^wS.(V>*»*i"V 


%      a  j 


f  ()>{)cQ(v  'H(J)     p(a)fn(v  )r  (J  Op'^a).  (IV-D-12) 

0  0       1  0     1     e°     1 

where  we  have  recognized  that 


Is  (r/)=L  (r\> 


^2'h(r^/)Lea(x2").  (IV-D-13) 


ea  *    ea 

This  last  expression  for  the  memory  operator  is  still  not  in 
a  particularly  convenient  form.  However,  there  is  a  transformation 
procedure,  outlined  in  Appendix  F,  that  will  cast  equation  IV-£-12 
in  a  form  which  we  can  more  easily  interpret.  This  technique 
puts  the  contributions  to  the  denominator,  which  are  in  the  form 
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of  integral  operators*  into  a  dynamically  shielded  electron-atom 
integration.     The  resulting  expression  is  IV-O-14 


M(co)  -_  U 
oj  A 

P 


2 


dr/dv/L0  (*  *)K»- ■* +iv  '•$  '-X     Ls  (?  -))""1 

i       i     ea     io)  i       i       ea  ea     i 

P 


P^^fo(v/)Lea(?/)p"1(a)'  (IV-45-14) 

where  the  dynamically  shielded  electron-atom  interaction  is  given 
by  equation  F-17: 


LD  (r  )= 
ea 


ea 


It  should  be  emphasized  that  no  approximation  was  made  in  going 
from  equation  IV-O-12  to  equation  IV-D-14.  The  operator  D~  (a,lc,Jc'ja>) 
includes  all  of  the  effects  of  the  operator  A  (to)  as  well  as  the 
effects  of  the  modified  Vlasov  operator;  its  functional  form 
may  be  inferred  from  equation  F-20.  Equations  IV-£-14  and  F-17 
give  the  expression  for  M(w)  correct  to  first  order  in  X  . 
YMle  equation  IV-D-14  is  exact  to  that  order  it  is  difficult  to 
evaluate  in  general.  However,  one  simple  approximation  to 
D""1(a,lc,,lc'jto)*  that  keeps  only  those  contributions  that  are 
zeroth.order  in  both  parameters  is  particularly  useful. 
Examination  of  equations  F-20  and  F-23  indicates  that,  in  this  approx- 
imation, D  (a,Jc,]c'ja))  becomes 

e(k;uvL  ) 
*   a 

E(kjto-L  )=[l-[dv,  k£-  f  (v  )C(k)].  (IV-D-15) 

J    or-L  -v  #k 
a  i 


which  gives  for  L  (r, ) 
ea 


LD  (x, )= 

ea  i 


D  /► 

L  (4) 


dke-ik'ri  -**- 


e(kja)-L  ) 
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eCkjw-I^)  will  be  recognized  as  a  generalized  dynamic  shielding 
function.  It  is  interesting  to  further  approximate  dynamic 
shielding  by  static  shielding,  giving  a  theory  for  which  all  of 
the  interactions  are  shielded  statically.  This  will  be  discussed 
further  in  the  next  chapter. 


CHAPTER  V 
DISCUSSION  OF  RESULTS 

V-A    Fully  Shielded  Unified  Theory 


In  Section  III-6  we  developed  an  expression  for  the  memory 
operator: 


M(  w)=~in 


d(l)L     (1)  7TT GLaa(l)p"1(a).  (III-B-28) 

ea        o)-LU;-V(a,l;a))       ea 


In  Section  IV-D,  V       and  V(a,l;w)  were  taken  to  first  order  in 

X       and  all  terms  in  M(co)  were  taken  to  zeroth  order  in  the 
ea 

plasma  parameter: 

X         r    .  L 


M(m) 

A     J      i     1     ea     l. 


=s-i  -Sa- 
to A 
P 


dr  dv  "L     (r  ')i(o>'-  -*■  +iv  '•$  '-1    r(i'  ) 
i     i     ea     l  u  i       i       ea  ea     i 

P 


♦inv/'f 'f.Cv/JCrUv  ')-Wrp(a)f(v')L    (f'Jp-Ha).         (IV-D-12) 
i        loi  o        i  o      i     ea      i 

We  then  proceeded  to  separate  the  operators  appearing  in  the 

denominator  into  two  parts:     those  containing  a  factor  of  the 

permutation  operator,  P. .,  and  those  not  containing  it.     The 

terms  containing  P.  .  were  combined  with  the  leading  L     (l)  to  form 
ij  ea 

a  dynamically  shielded  interaction   (Appendix  F).     The  remaining 
terms  containing  no  factor  of  P. .  were  left  in  the  denominator  to 
effectively  shield  the  electron-atom  interaction  appearing  there: 

w  p  J  "p 

P(«)f,(v1')lJa(?  Op"^.).  (IV-D-14) 


TD     /■*■     \- 

ea      i 


dfdl'L    («^)D"1(a,1c^';a))e"i,r*?1 


ea 
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where  the  operator,  D  (a,k,k  i^),   contains  all  the  effects 
included  in  A(w).  In  this  section  we  discuss  some  of  the 
consequences  of  the  functional  form  of  equation  IV-45-14. 

This  result  is  interesting  because  it  leads  to  a  unified  theory 
in  which  all  of  the  interactions  are  shielded.  The  unified  theory 
of  VCS  did  not  explicitly  include  shielding,  while  the  unified 
theory  of  Capes  and  Voslamber  led  to  an  expression  for  M(w)  in  which 
only  those  interactions  appearing  in  the  numerator  were  shielded. 
The  weak  coupling  limit  used  to  develop  equation  IV-£-l2  is  more 
systematic  than  that  used  by  either  VCS  of  Capes  and  Voslamber  and 
leads  to  the  fully  shielded  result  presented  here. 

As  indicated  in  Section  IV-©,  a  further  approximation  would 

be  to  take  the  shielding  function  D~  (a,k,k"so))  to  zeroth  order 

on  both  A  and  X     .  This  leads  to  a  result  similar  to  that  of 
ea 

Capes  and  Voslamber  (Section  II-C),  except  for  the  static 
shielding  appearing  in  the  denominator.  As  a  further  approximation 
to  this  last  result  we  replace  dynamic  shielding  by  static 
shielding.  This  yields  an  expression  for  the  memory  operator 
in  which  all  of  the  interactions  are  statically  shielded  in  the 
same  way: 


M(w)~in 


,s  /-*■ 


d(i)Lja(x/)i(«an;1^1-tja(jx))-4p(«)f0(v1)L;aorxOp"i(«). 

■ 
Lea(^)=Lea(?i)^Jd(2)fo(v2)h(?^)Lea(?2)'  (»-**) 

where  h(r  r  )  is  the  electron -electron  pair  correlation  function. 
The  fact  that  all  interactions  in  equation  IV-A-1  are  the  same  allows 
us  to  integrate  by  parts,  following  the  procedure  of  VCS  outlined 


in  Section  II-*.     Applying  this  procedure  yields 


f  <-    + 


M(w)=-in(w-L  )  dTjdv  [IKO-Uf^v,  )(w-L ), 


U^^exp{-; 


dt'IT  (r)}*  (V/^A-2) 


0        ea 


where  U(t)  is  the  time  development  operator  in  the  interaction 
representation  and  T  is  the  time  ordering  operator.     Equation 
IV-A-3  closely  resembles  the  expression  which  VCS  evaluated,  except 
that  now  the  electron-atom  interaction  appearing  in  the  exponential 
is  shielded.     VCS  evaluated  the  integral  in  equation  IV-A-2  by 
neglecting  time  ordering,  assuming  that  electrons  followed 
straight  line  trajectories,  and  cutting  off  the  spatial 
integrations  at  the  Debye  sphere.    The  numerical  results  which 
they  obtained  using  their  method  have,  in  general,  agreed  well  with 
experiment  in  spite  of  the  ad  hoc  way  in  which  they  accounted  for 
electron  correlations.     Since  an  integral  over  a  statically 
shielded  interaction  can  frequently  be  approximated  with  good 
accuracy,  by  a  bare  interaction  and  some  sort  of  cutoff,  we 
can  understand  why  the  procedure  employed  by  VCS  was  so  successful. 

V-B     Conclusion 

In  Chapter  III  we  developed  a  general  expression  for  the 
line  shape  function  for  a  dipole  radiator  immersed  in  a  one 
component  perturbing  fluid.     The  calculation  was  fully  quantum 
mechanical  and  it  was  never  necessary  for  us  to  specify  the  nature 
of  either  the  perturbing  fluid  or  the  radiator.     A  formally  exact 
expression  for  the  line  shape  was  obtained  in  a  form  suitable  for 
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approximation.  Then  in  Chapter  IV  we  assumed  classical  perturbers 
interacting  with  the  radiator  via  a  Coulomb  potential,  and  we 
applied  the  weak  coupling  limit.  It  was  determined  that  the  weak 
coupling  limit  could  be  applied  in  a  variety  of  ways  and  we 
discovered  that  we  could  reproduce  several  of  the  existing  theories 
of  line  broadening,  depending  on  how  these  limits  were  used.  We 
also  developed  a  unified  theory,  which  was  fully  shielded,  by 
expanding  M(w)  to  zeroth  order  in  the  plasma  parameter. 

We  should  emphasize  the  generality  of  the  method  employed 

35 
here.  While  we  have  used  some  of  the  ideas  of  Gross  *  in  our 

development  of  the  short  time  limit,  this  work  was  primarily 

32  33 
based  on  a  technique  developed  by  Mazenko  *   to  study  the  general 

problem  of  time  correlation  functions.  The  development  of  Chapter  III 

was  similar  to  a  study  of  the  velocity  correlation  function  made 

32 
by  Mazenko   except  that  we  have  two  different  kinds  of  particles: 

perturbers  and  radiators;  hence  we  have  two  different  kinds  of 

interactions.  In  Chapter  IV  we  restricted  our  consideration  to 

a  plasma  for  which  the  classical  path  approximation  was  valid, 

and  applied  a  weak  coupling  limit  based  on  the  nature  of  these 

interactions. 

We  should  note,  however,  that  it  is  ,  in  general,  possible 

to  parallel  more  closely  the  work  of  Gross  or  Mazenko  and  bypass 

the  weak  coupling  limit.  In  Section  III-C  we  discussed  the  use  of 

short  time  limits  to  formally  close  the  BBGKY  hierarchy  of  kinetic 

equations.  Hopefully  it  will  be  possible  using  this  approximation 

procedure,  to  evaluate  a  fully  quantum  mechanical  expression  for 

the  line  shape  without  applying  the  weak  coupling  limit.  Such 

an  expression  would  be  necessary  to  properly  treat  very  dense 

45 


plasmas,  such  as  laser  produced  plasmas, 
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APPENDIX  A 
EQUATION  FCR  K(a,lj  ) 

In  this  appendix  the  second  equation  of  the  hierarchy,  equation 
IIIiA-10,  will  be  put  in  a  form  that  allows  us  to  formally  solve 
for  the  operator,  K(a,l;co),  and  in  turn,  to  arrive  at  a  formal 
expression  for  the  memory  operator,  M(to).  We  have  used  equation 
I II»© -12,  which  separates  5(a,l;t)  and  D(a,l,2jt)  into  their 
short  time  functionals  of  D(ajt)  and  remainder  terms,  p(a,ljt)  and 
?(a,l,2jt),  respectively.  The  result  obtained  by  substituting 
this  separated  form  into  III-^A-10  is 

(A-  +iL(l))P(a,l)p-1(a)DW;)+  (^  +iL(l))?(a,l;t) 

=-inTr  {(L  (2)+L2p1)p(a,l,2)p""1(a)}D^(a;t) 

-4nTr2{(Lea(2)+Lee)P"(a,l,2jt) }.  (III-B-22) 

We  may  simplify  equation  III-©~22  by  recalling  the  first  equation 
of  the  hierarchy 

(■&£•    +iLa)rT(a;t)=-inTri{Lea(l)Dv(a,lSt)}.  (III-A-IO) 

With  the  separated  expression  for  Dv(a,ljt)  this  equation  becomes 

■5rU(ajt)=^.LD(a;t)-inTx  {L  a(l)p(a;t)p"'1(a))Dv(ajt) 

a  u  a  i  e  a 

-inTr  {L  (l)?(a;ljt)}  (A-l) 
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which,  when  substituted  into  equation  III-43-22,  yields 
(3_+Lf>))P(a,l;t)nnTr.{(L  ,(2)+L")P(«,l,2jt)> 

j\  x  <-    6  3       66 

-ip(a,l)p"*1(a)nTr.{L  (l)P(a,l;t)}=^-i  [L(l)p(a,D 

i   Q  3 

^Tr2{(Lea(2)+Le^:f<a,l,2)}]p"1(a)D(ajt) 
+ip(a,l)p"1(a)  [Lap(a)-^Tr1{Lea(l)p(a.D  }]p"1(a)D(ajt).   (A-2) 


The  two  terms  in  brackets  on  the  right  hand  side  of  equation  A-2 
can  also  be  simplified  by  looking  at  the  first  two  equations  of  the 
equilibrium  hierarchy I 

[I  p(a)«nTx.{L  B(l)p(a,l)>]p"1(a)5(a|t) 

^l...N{LPaN}p"1(a)iJ(ait)  (A"3) 

[L(l)p(a,l)^Tr  I(L  (2)-Lfo1)p(a.l,2)}]p"1(a)^(ajt) 

/.         63       66 

=TV..N{LPaN}p'"1(a^(ajt)-  (A-4) 

It  will  be  noted  that  Lp  «  would  vanish  except  for  the  fact  that 


it  operates  on  a  function  of  atomic  coordinates,  p  (a)D(a;t). 
The  operators,  L  and  p^,  will,  however,  commute;  that  is, 
Ip^,^  NL,  and  we  rewrite  equations  A-3  and  A-4, 

[L  pCaJ^Tr.R  (l)p(a,l)}]p""1(a)^(a;t) 

3  1   63 

=  [p(a)L  •»nTr1{p(a,l)L  (l)}]P"1(a)^(aSt)  (A-5) 

3      1  6  3 

[L(l)p(a,l)-»nTr o{(L  a(2)+L")P(a,l';2)}]p"1(a)U(a;t) 

2    63       66 
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=  [p(a,l)(La+Lea(l))+nTr2{p(a,l,2)Lea(2)}]p"1(a)^(ajt)     (A-6) 

where  we  have  also  made  use  of  the  fact  that  terms  with  L  (j)  and 
L  ,  operating  on  purely  atomic  operators,  vanish.  If  equations 

66 

A-5  and  A-6  are  substituted  into  equation  A-2  we  get 
(it  +iL(l))^a*ljt)+inTr2{(Lea(2)+Lee)p(a,l,2;t)  } 
-inpCa^Dp"^)^^!  (l)p(a,ljt)}=-ip(a,l)L  (Dp-^ajDCajt) 

1   6  9  ©3 

-inTro{p(a,l,2)L  (2)}p~1(a)^(a;t)+inp(a,l)p"'1(a)x 
xTr{p(a,l)L  (Dlp^Ca^Cajt).  (A-7) 

1  69 

Equation  A-7  is  useful  because  with  it  we  have  an  equation 
for  f(a#l;t)  and  P"(a,l,2jt)  in  terms  of  JJ(ast).  We  can  define 
a  time  independent  atomic  operator,  GL  (l),  given  by 


ea 


GL     (l)p   '(a^astJ^Ca^l)!     (l)p**1(a)D(ajt)-^Tr9{p(a,l,2)L     (2) 

xp~Ha)tf(ajt)^p(a,l)p~Ha)Tr>(a,l)L     (lnp^aMajt)       (A-8) 

1  69 

with  which  we  can  rewrite  equation  A-7  in  a  more  compact  manner; 
(■^  +iL(l))P(a,ljt)+inTr2{(Lea(2)+L^)?(a,l,2;t) } 

-inp(a>l)p""1(a)Trn{L    (l)p(a,ljt)}=-iGL     (iJp'^ajDCajt).     (A-9) 

16  3  6  9 

It  is  now  possible  to  Laplace  transform  equation  A-9  with  the  result 
(ar.L(lV(a,lju>)-nTr   {(L     (2)+L     )pW,2;  J  } 

2  6  3  66 
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^np(a,l)p     (a)Tr,{L    (l)p(a,l5u)}=GL  ^(l)p"  (a)D(a;to).  (A-10) 

16  3  63 

From  equation  III-B-16  we  recall  that  P(a,l*,*sjw)=K(a,l,,,sjw)0(aja)) 
hence  we  are  able  to  get  an  equation  for  K(a,ljw); 

(o)-L(l))k(a,l;o))-inTr  {(L     (2)+L^)R(a»l,2ju))} 

-npUjlJp^UJTr.a     (l)K(aal;a))}=GL     (Dp^a).  (A-ll) 

j.     ea  ea 

Finally  in  this  appendix  we  will  examine  the  operator  GL    (1), 

defined  by  equation  A-8,  noting  that  all  of  the  terms  oh  the  right 

hand  side  of  that  equation  contain  a  factor  of  the  electron-atom 

interaction,  L     (j).     If  we  define  a  permutation  operator.  P.., 
ea  xj 

that  will  change  the  functional  dependence  of  all  operators 
appearing  to  its  right  from  functions  of  particle  j  coordinates  to 
functions  of  particle  i  coordinates,  then  we  can  rewrite  equation  A-8: 

<*    pJ(a)3(ajt)Ma,l)L     (lJp^Ca)     ^Jt) 

^Tr?{p(a,l,2)P21)L     (l)p~  (a)D(ajt) 

-np(a,l)p"1(a)Tr2{p(a,2)P21}Lea(l)ip-\a)^(a;t).  <A-12) 

Thus  we  are  able  to  define 

G=p(a,l)+nTr2{P(a,l,2)P2i}-p(a,l)p"1(a)nTr2{p(a,2)P2i}  . 

(A-13) 


APPENDIX  B 
ANALYSIS  OF  V(a,ljoi) 


The  purpose  of  this  appendix  is  to  validate  equation  III-C-26, 

V(a,ljjK(a,lja))=nTr  {L  (2)+L  )R(a,l,2ju)) 

2  ea    ee 

-pUDp-HaVTr  {L     (l)K(a,lj<J},  (lII-B-26) 

1    ®  *• 

and  to  study  the  short  time  (infinite  frequency)  limit  of  V(a,ljco). 
Our  first  step  is  to  express  K(a*l>2jaj)  as  some  functional  of 
K(a,l;o)).  Following  Mazenko  *   we  will  use  exactly  the  same 
technique  for  expressing  K(a,l,2jco)  as  a  functional  of  K(a*l;w) 
that  we  used  to  express  l3(a,ljt)  as  a  functional  of  ^(ajt). 
Equations  III-*-5  and  III-3-6  suggest  that  we  should  attempt 
to  express  K(a,l,,#s;u))  as  a  linear  map  of  a  function  spanning  the 
space  of  the  coordinates  of  the  atom  and  particle  1  onto  the 
space  of  the  coordinates  of  the  atom  and  s  particles. 

First  we  recall  the  definition  of  K(a,l,,#sju)),  equation 

ni-B-n; 

K(a,l^»sja))=U(a,l"*sjco)0"1(aja))-p(a,l»*'s)p"1(a).   (III-B-17) 

From  equations  HI-£-l  and  III-B-3  we  have  the  definitions  of 
U(a,l-»»sjco)  and  0(aju): 

0(a,l— siuKr^^..^^  ^  PaN}p"1(a)  (HI-B-l) 

U(aj  J=Tr1###N{/  ^  paN}p"1(a).  (III-B-3) 

Making  use  of  the  fact  that  p  „  and  L  commute,  and  substituting 
equations  IIIH3-1  and  III-e-3  into  equation  III-B-17  we  arrive  at 
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K(a,l...SjJ=Txfi+1..#N{VNp  m  £-  }p~l  (aJiTr^.^/  ^  P  ^p"1  (a)]"1- 

-p(a,l«»«*)p"1(a).  (B-l) 

Noting  that 

^(a'w)=  0,-L  hU(o,)  p(a)J  (IIM-19) 

a 

a 
We  can  rewrite  equation  B-l  in  a  more  useful  form: 

K(a,l...s;a))=-iTrs+1...N{PaN  ££■  }p"1(a) [aHLa«e-M(M) ]-o(a.l—«)p"1(a)-  (B-3) 

where  B  and  M(w)  are  given  by  equations  III-B-19  and  III-B-28 
respectively.  We  now  make  use  of  the  identity, 


u« 


1-  =  -! — _L.i(I^L)-*  (b^i) 

)-L   w-L    co-L     a  w-L 
a  a 

in  order  to  get 

K^-l— ^)=-iTrs+1...N{paN[l-^L  i  zL„(J)]pT  x 

J         a 

xp"1(a)[co-L  -B-4l(a))]-p(a,l*'«s)p",1(a)j  (B-5) 

a 

21 

we  have  observed  that  the  result  of  L •  (j)  and  L   operating  on  purely 

e  a       e© 

atomic  operators  is  zero.  It  would  be  tempting  to  cancel 


-^r-x(co-L  ),  but  we  cannot  since  L  will  not  in  general  commute 
a-L     a  a 

a 

with  p(a)=Tr1#>#N{p  „}.  However,  from  the  first  equation  of 
the  equilibrium  hierarchy,  equation  A-5,  we  know  that 

PJ    (a)[L  •^]=p"1(a)[Lp(a)-*nTr1{L  ADpUmDUp^U) 

3  3  1        6q 

=p""1(a)[p(a)L  +nTr{p(a,l)L     (l)}]p"1(a) 
a         l  ea 

=Lp"1(a)^",1(a)nTr  {p(a,l)L     (DJJp^ta).  (B-6) 

a  i  ea 
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If  this  is  substituted  into  equation  B-6  we  get 

K(a,l-s;-)=.iTVl#..N{paN^rZLea(j)}p-1(a) 

+iTVl...M{p*  ^T>p"1(a)[nTr1{p(a,l)Lea(l)p"1(a)+M(u))].     (B-7) 
We  now  consider  equation  III-B-28 

M(a>)=-inTr  {L     (l)  rrr-4 GL  a(l)p~1(a), 

1     ea        w-LU,-VUl;u>)       ea 

which  we  will  substitute  into  equation  B-7  and  which  is  consistent 

with  the  formal  definition  of  V(a,l;w)  in  equation  III«*-26.  Now 

we  can  extract  a  factor  of  L  (l)p""  (a)  from  all  of  the  terms 

ea 

in  equation  B-7  by  introducing  the  permutation  operator  defined  in 
Appendix  A.     Hence  we  are  able  to  define  the  operator,  X(a*l,,#sju>). 

K(a,l«*»Siu))=X(a/il,,,S|u)L     (l)p"  (a) 

X(a,l---sju)Lea(l)p"1(a)=-iTVl...N{paN^rEPji}Lea(l)p--1(a) 

wvi..*tp*  ^r  '>p",1(«)p»Tr2{p(«.i)p21}ibt(i)p"1u) 

-inTr  {L    (l)  rry-i GL     (l)p"1(a)}.  (B-9) 

2    ea         arLU)^(a,l*tt)         ea 

We  will  now  assume  that  X(a,lj w)  has  an  inverse  and  write 

K(a,l»«»sju,|K(a,l;a)))=X(a,l»».sja))^."'1(a,l;(Jj)K(a,ljJ  (B-10) 

which,  when  substituted  into  equation  IIIHB-26,  yields 

V(a,l;u)K(a,l;u)=*iTr  {(L  (2)+L ■)XU»Uu)^Tl(a»lia)  }K(a,l;J 

-pU^lJp^UJnTr  (Lea(2)P  }R(a,ljJ 

V(a,ljJ=%iTr  {(L  (2)+L ■  )^(a,l,2;  J^"!(a,l;  J  } 

0     6d       cc 
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-p(a,l)p"1(a)nTr2{Lea(2)P21}.  (B-ll) 

Thus,  we  have  justified  equation  III-B-26  and  in  doing  so  we  have 

derived  an  exact  formal  expression  for  V(a,l;w). 

We  will  now  separate  V(a,lj^)  into  a  frequency  dependent  and 

independent  parts 

V(a,lju>)^(a,l;o*=~)-HV  (a,ljw);  (B-12) 

c 

V(a,lj^»)=hTr2{(Lea(2)+Le2e1)?(:(a,l,2;a^«>)^"1(a>ljaf^)} 

-*>(a,l)p"  (a)nTr?{L  (2)P  },  (B-13) 

Vc(a,ljco)=^Tr2{(Lea(2)+L^)5ic(a,l,2;aJ)?C"1(a,lja)) 

-fc(a,1.2ju*»)X"*(a,lso)=*»)}.  (B-14) 

It  is  possible  to  further  analyze  the  expression  for  VCa,!;^30) 
by  realizing  that  £im  ~t*  =  *■  •  Neglecting  terms  that  are 
higher  than  first  order  in  *•  from  equation  B-9  allows  us  to 
write 

X(a,l-..s;W-)Lea(l)P^(a)=^[Trs+1...N{paNEPjl} 

ip(a,l—*)p"l(a)Tr,{p(a,2)Poi}]L  (Dp-^a).       (B-15) 

When  this  is  substituted  into  equation  B-10  we  get 

X(a,l--Sja)=»)^l"J(ajl;a)=^)K(a,lja))=[Trs+1_#N{paN  ePj   } 

-np(a,l-»s)p"1(a)Tr2{P(a>2)P21}][Tr2.#.N{PaN  jP^} 

-np(a,l)P   1(a)Tx2{P(a.2)P21}]'1k(a.l|(o),  (B-16) 

where  it  will  be  noted  that  the  permutation  operator,  P. .,  operates 
on  everything  that  appears  to  its  right.  Hence,  we  are  able  to 
write  V(a,ljaia»)  explicitly  in  terms  of  known  functions,  while 


75 


the  formally  exact  collisional  term,  V  (a^l;^),  requires  the 

c 

next  equation  of  the  hierarchy. 


APPENDIX  C 
KINETIC  EQUATION  FCR  x(a,l,2;aj) 

In  this  .appendix  we  will  transform  the  third  equation  of  the 
hierarchy, 

(  |-  +iL(2))D(a,l,2jt)=-inTr3{(Lea(3)+L^+Le32)D(a,l,2,3jt),  (III-A-10) 

into  an  equation  for  x(a*l*2$w)  defined  by  equation  III-  B-36.  As 
in  Appendix  B  we  use  equation  III-8-12  to  get 

(  |^+iL(2))p(a,l,2)p-1(a)D(a,t)+(  ~  +iL(2))P(a,l,2;t) 

=inTr,*(L    (8)+L    +L    HpiaA&Sh"1  U)Ua}t) 

3       ea  ee     ee 

+P(a,l,2,3;t)].  (Ol) 

From  the  first  equation  of  the  hierarchy  one  obtains 

~r^(ait)=-iLr5(ajt)-inTr  {L     (l)p(a,l)}p"1(a)B(a;t) 
dx  a  i     ea 

-inTr  {L     (l)^(a,ljt)},  (III-A-10) 

1     ea 

which  can  be  substituted  into  equation  C-l  to  give 
(  2-  +iL(2))?(a,l,2it)-inTri{Lea(l)?(a,l;t)} 


+inTr  {(L    (3)+Lfi+L")?(a*l,2,3jt)=^i J2^p(a,l,2)p'"1(a)D(ajt) 

3       ea  ee     ee  r  r 

+p(a,l,2)p"1(a)Larj(ait)+ip(a>l,2)p-1(a)nTri{Lea(l)p(a,l)} 
P**1(a)D(ajt)-inTr,{(L     (3)+L"+L3i)p(a,l,2,3)}p"1(a)U(a;t). 


3lN  ea  ee     ee 

(C-2) 
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The  Laplace  transform  of  equation  C-2  yields 
(orL^jKUl^jJ+nTr.tL  a(l)K(a,lj<J 

HiTr,{(L    (3)+L"+L")K(a,lA3ja)))=L^)p(a.l»2)p-1(a) 

-p(a,l,2)p-1(a)La+nTr3{(Lea(3)+Le3e1+Le3e2)x 

xp(a,l,2.3)}p-1(a)-p(a,l,2)p"1(a)nTrl{L ^^  (l)p(a,l)}p"*1(a),     (C-3) 

where  we  have  explicitly  used  the  fact  that  PCa^l'^'sj^O)^  and 
the  relation  P(a,l,,,sja3)=K(a,l*,*Sjoo)D0aja))*  equations  IIIH3-13 
and  III-&-16,  respectively.     It  is  now  useful  to  take  the  inverse 
Laplace  transform  of  equation  C-3* 

(  5~+iL(2))K(a,l,2jt)-ip(a,l,2)p'"1(a)nTr3{Lea(3)K(a,3;t)} 

=-inTr3{(Lea(3)+Le3^Le3e2)K(a,l,23;t)},  (C-4) 

where  we  have  noted  that 

K(a,l,2;t=0)=5-p(a,1.2)p"1(a)iLa+iL(2)p(a,1.2)p"1(a) 

nnTr3{(Lea(3)+L^H-Le3|)p(a,l,33)}p-1(a) 

-in(a,l,2)p-1(a)nTr1{Lo  a(l)p(a,l)}p"1(a).  (OB) 

Equation  C-4  gives  a  relation  between  K(a>l;t),  K(a,l,2jt)ji 
and  K(a*l*2,3;t)  which  will  be  closed  in  a  manner  analogous  to 
that  shown  in  equation  III-e-^6.  We  start  by  using  equation  III-£-33 

K(a,l»»»s;t}K(a,ljt))=K(a,l*'*sjt=0|K(a,ljt))+X(a,l«"s;t)     (C-6) 
where 
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K(a,l"*s;t=0|K(aJljt))=^(a,l'"s;ar^)?C"1(a,lja))K(a,ljt)  (C-7) 

Again  the  correspondence  between  the  t+0  limit  and  the  cj-*»  limit 
should  be  noted.  Equation  C-7  may  be  substituted  into  equation  C-4 
with  the  result  that 

(  ^•+iL(2))^(a,l,2jo^)?C""1(aJlja)  =°)K(a,l;t) 

+(  |^  +iL(2))x(a,l,2it)-p(a,l,2)p-1(a)nTr3{Lea(3)K(a,3;t)} 

— inTr3{(L  (3)+L"+L")  [^(a,l,2,3jaF=°°)X"1(a,l;^=-)x 
3   ea    ee  ee 

xK(a,l;t)+X(a,l,2,3jt)].  (C-8) 

Using  the  inverse  Laplace  transform  of  equation  A-ll,  we  now 
eliminate  rrK(a,ljt)  from  equation  08: 

(  |-  +iL(2))X(a,l,2jt)HnTr3{(Lea(3)+Le3e1+Le3e2)X(a,l,2,3St)} 

-iX(a,l,2jaR»)X"1(a,lsaF~)nTr,{(L  a(2)+L21)x(a,l,2jt)} 
=i-iL(2V(a,l,2;a^)?<."1(a>liaF-)K(a,lst)+iX(a<l,2;aF»)x 

xX"1(a,l}aP»)L(l)K(a#lit)+p(a,l,2)p"1(a)nTr,{L  (3)K(a,2;t)} 

-Ma,l*2soP»)^"1(a,l;aF-)p(a<l)p"1(a)nTr3{L  (3)K(a,3;t)} 

g  a 

+inTr  J(L  (3)+L3i+L3o2)[%(a,l,2-3;aF«-)^"1(a,ljaF»)K(a,ljt)} 

nX(a*l,2;aF-)^"1(a,l;oi=oo)nTr2{(Lea(2)+Le2e1)X(a,l,2;aF«) 

X~1(a,ljar»)K(a,ljt)}.  (0-9) 

Since  X(a,l»»»sj  w)=x(a,l**,s;  w)K(a,lj  w),  the  Laplace  transform  of 
the  above -.equation  is 
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(arL(2))x(a»l,2jco)+^(a,l,2iaP»)^"1(a,l;ar<»)x 
xhTx2{(Lea(2)+L^)x(a,l,2iW)HiTrs{(Lea(3)+L»+L^)x 

xx(a,l,2,3joo)}=S(a,l,2)  (O10) 

where  S(a«lA2)  is  a  frequency  independent,  inhomogeneous  source  term 
given  by 

SCa^l^jKCa^ljwJ^XCa^^aj^^-^a^liaf^JfL^^Ca,!^) 

-p(a,l)p"1(a)nTr2{La(2)K(a,ljw)}-hf»Tr,{(L  a(2)+L")x 

x^(a,l,2Sa)^)^:"1ta,l;a3=»)K(a,ljw)]+L(2W(a,l,2ja3=b0)x 

x^"*1  Uljw^KUl^Majl^p"1  (a)nTr,{L  (3)K(a,3;co)> 

9        6  3 

%Tr3{(Lea(3)+Lee+Lee)^(a,1.2,3ja>»o)X""1(a,liu!*>)x 
xK(a,ljw).  (Oil) 

This  is  the  result  needed  in  order  to  solve  for  V  (a,ljw). 

c 


APPENDIX  D 
DIMENSIOIIESS  LNITS 


In  order  to  understand  the  justification  for  a  plasxaa 
parameter  expansion*  let  us  first  look  at  the  linearized  Vlasov 
equation  for  a  one-component  plasma: 

<  ft  ^r^)fU)^Vt)=m  V*,Cvx)fd(2)  -^f(l)(r2,v2St),(D-D 


3vi 


*i 


where  <f>2f  -*—  is  the  interparticle  potential.  In  the  literature 

it  has  been  shown  that  the  natural  units  should  be 

i  i. 

w  •~(4iine2A)2j     L^KI^Tine2)2; 


u^o-bflL,     x->x'=f-.     v^v-b-JL..  (D-2) 

Wp  *D  *Dup 

If  equation  D-l  is  scaled  with  these  dimensionless  units  the  result 
is 


({p^I-^)«(i)(?1.*1.t)-«jij^-yo)(vl-) 

P     3Vj 


-*■  ^    r 
9r  12 


x  f(l)(?  'v  '.t)=i-^--f  (v  ')fd(2')  -3 MX)(r  ,v  St).(EK3) 


"I        2     "  4TT    ^        o 


3r   '    r' 

1  12 


There  are  two  observations  that  are  worth  making  about  equation  D-3. 

The  first  of  these  is  the  fact  that  whenever  L2*   appears  in  a 

kinetic  equation  which  is  scaled  to  the  dimensionless  units 

defined  by  equation  D-2,  it  will  usually  be  accompanied  by  a  factor 

of  -*- •     The  second  observation  is  that,  even  though  the  right 
w 
P 
hand  side  of  the  Vlasov  equation,  equation  D-l,  is  first  order  in 

the  electron-electron  coupling  strength,  e  ,  the  right  hand  side 
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of  equation  D-3,  when  expressed  in  dimension  less  units,  is 

zeroth  order  in  this  parameter.     Since  this  term  will  not  be  small 

for  long  range  forces,  this  suggests  that  perhaps  e2  is  not 

the  natural  expansion  parameter.     Let  us  examine  the  following 

expression: 

up    ee       una    rx  v,    r2    v2  r12        ^        r,    Vj    r2    v2    r^ 

4™V    ?!       *,       ?2      ^2ri2  Xl       Vl      r2      V2X 


12 


fc_l 


where  we  identify  A=  V^y J  as  the  plasma  parameter.  For  typical 
experimental  temperatures  and  densities  this  parameter  is  about 
0.05;  this  suggests  that  it  might  be  a  good  expansion  parameter. 
It  should  be  noticed  that,  as  in  the  case  of  the  linearized 

Vlasov  equation,  the  factor  of  *j-  L   will  often  appear  under  a 

P 

three  dimensional  spatial  integral  and  be  multiplied  by  a 

factor  of  the  density.     In  this  event,  the  scaling  to  dimension  less 
units  produces  a  term  which  is  zeroth  order  in  the  plasma  parameter 
but  first  or  higher  order  in  the  electron-electron  coupling  strength; 
thus,  in  general,  a  systematic  expansion  in  the  plasma  parameter 
will  include  different  terms  than  a  systematic  expansion  in 
the  electron-electron  coupling  strength.     The  advantages  of  a 
plasma  parameter  expansion  over  an  expansion  in  the  coupling 
constant  for  systems  with  long  range  forces     are  discussed  in 
most  texts  on  plasma  physics.     * 

It  Will  now  be  useful  to  rewrite  some  of  the  quantities  that 
appear  in  the  line  broadening  formalism  of  Chapter  III  as  functions 
Of  the  dimensionless  variables  defined  by  equation  D-2: 


oz 


^a*w)=u>-LHUu)    Pa^ 


w  D(ajco)- 


l< 


p(a)d\ 


(lII-G-17) 
(D-5) 


p  ,   1    L     MLwL 

Note  that  both  B  and  M(^)  appear  with  a  factor  of  "JJ" 

P 

d(l)p(a,l)  -" p    (a) 

P 

L     (1)                t  UD     »i 

jaLi rk G-" — P    (a) 


CO 


R* 


d(l) 


CO. 


b) 


"TT5  -  w 

* y,(q.lj")  p 


(D-6) 
(D-7) 


-p        -p 

where  the  trace  operation  of  equations  IIIH3-25  and  III-e-26  has 

been  replaced  by  integrals  indicating  that  the  classical  path 

approximation  has  been  made.     Again,  note  that  every  factor  of  L     (j) 

ea 

appears  with  a  factor  of  jj- which  must  be  taken  into  account 


when  this  operator  is  scaled 
L  (1)     2   r 


-§£. 


to, 


2 

e 
P 


?        e        R 

•R=  — ?-r    r-  • 


r. 


IS. |»    "Vd    *     |J,|«! 


(D-8) 


If  we  take  R  to  be  approximately  the  Bohr  radius,   a  ,  then  we  can 
define  the  electron-atom  coupling  strength  to  be 


X     = 
ea 


u 


e2a. 


& 


(D-9) 


which,  for  typical  experimental  conditions,  is  about  0.01;  hence, 
we  can  expect  it  to  be  a  good  expansion  parameter. 

It  is  now  useful  to  examine  G  which  is  given  classically  by 


L     (1) 


L     (1) 


«f PJ  (a)=p(a,l)  -f p'HaW  d(2)p(a,l,2)  -**~  p"  (a) 


10 

P 
-np(a,l)p"  (a) 


L     <2>     - 


L     (1) 


*p 


d(l)p(a,l)  -" 


-l 


w_ 


(a). 


(D-10) 
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The  simplest  way  to  deal  with  p(a,l#"s)  is  to  define 
p(a,l,,'s)=p(a)p(l),,'p(s)g(a,l"*s) 
p(j)rf0(Vj)A 


30 


(D-ll) 


where  f  (v.)  is  the  Maxwell-eoltzman  velocity  distribution.  We  will 
o  J 

not  scale  g(a,l***s)  here;  but  it  should  be  noted  that  to  zeroth 
order  in  \    ,  ${a»Vs)-*q{V$)»   where  g(l»««s)  is  the  usual 

spatial  reduced  distribution  function  which  may  easily  be  expanded 

44  /   \ 

in  the  plasma  parameter,  A.         Since  L     \j)  is  a  function  of 

r.  only,  equation  D-10  becomes 


1  JSi')     ., 


J& 


f  (viO 


w. 


p"(»)=  ^r-  p(«)[»(-D>,A.(S1')+  T 


dr  'x 

i 


x(gUl,2)-gUl)g(a,2))Lea(r2')] 

and  when  this  is  substituted  into  equation  D-7,  that  equation 
becomes,  in  dimensionless  coordinates, 


(D-12) 


ml  -. 


w. 


to. 


dv'dr  'L    (!') 
i      i     ea     i 


O) 


—J : f  (v  ')p(a)x 


w. 


0). 


x[g(a,l)Lea(r/)+-^-  fdr2(g(a.l.2)-g(a,l)g(a,2))Lea(^p"ll(a)XD-13) 

All  that  remains  now  is  to  express  1/   '+V(a,l;w)  in  dimensionless 
coordinates.     Using  equation  III-B-8  we  get 


(1) 


iv 


$        L    (1) 


jl«  ~j — j-  +  -aa 


03 


03 


03 


0) 


=  -a-~ivv--vr'+JL  L     (x    ')• 
oj  i     i    ea  ea     i 


(D-14) 


p  P  P  P 

The  operator  V(a,l;  <*)  will  be  treated  in  two  parts:   the  mean  field 

part,  vCa*!;"3")*  and  the  collisional  part,  V  (a,l;oj).     The 

former  was  given  by  equation  III-B-30 

V(a,li^)=nfd(2)(L     (2)+L21)[p(a,l,2)(l+P21) 
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+n 


d(3)p(a,l,2,3)P„  -np(a,l,2)P   Xa) 


31 


d(2)p(a,2)P     ]X 

21 


x[p(a,l)-*n 


-1/ 


d(2)p(a,l,2)P3l-np(a,l)p     (a) 


d(2)p(a,2)P2il' 


rl 


-p(a,l)p"\a)  d(l)L    (l)P     . 
I  ea         2i 


(D-15) 


Using  equation  D-ll,  the  above  may  be  scaled  to  give 


y,(,g  J;.<r.°°}  s 


^x'dv   ^Leatf2')4i(Ar:;ir  ).  ^){p(a)f  (v   -)* 0(v2')x 

p  J  <ir  12  3V 

:[g(a,l,2)(l+P21)+Wd(3')f0(vJ(g(a,lA3Wa,l,2)g(a,3))P3l]}x 


x{p(a)f0(v1')[gUl)+*-   d(2')f  (v    1(g(a,l,2)-g(a,l)g(a,2))P21]}' 


rl 


-Ag(a,i) 


A  ' 


^W'v 


(D-16) 


V  (a>l|to)  defined  by  equation  IIIH3-«39ji 
c 


V  (a,lju>)=i-in 
c 


d(2)(L     (2)+L     ) 


uhl*    "4na,l*2jco; 


S(a,l,2)     (III-6-34) 


is  more  difficult  to  scale.  It  can  be  shown  that  all  terms  in 

S(a,l,2),  given  by  equation  C-12*  contain  a  factor  of  L  M   L  , 

ea   ee 

or  (g(l2)-l)  which  does  not  appear  under  an  integral?  thus  this 

term  is  always  first  order  in  either  A   or  A.  In  dimensionless 

ea 

coordinates,  equation  III-B-39  becomes 

V  (a,l;ea)         t  A  a         i  ^ 

(2)         "  l 

i(B'-Xi EUi*2wl)-is(a,i,2)i 

(0  0) 

P  P 

T(2)       L 

i s-A^S  -.$  '-iv"  --$  ^^4+,^.^h     L     (r  ')+A     L     (?  '  ) 

a  u  1       1        2       2  rj     wl     r^     v2       ea  eav  j         ea  ea    2 

(D-17) 


showing  that  V  (a,l;w)  is  at  least  first  order  in  Aea  or  A.  Since 
no  real  attempt  is  made  to  study  Vc(a,ljto)  in  this  dissertation, 
other  than  to  identify  its  order  in  an  expansion  in  Xeg  and  A,  we 
will  not  attempt  to  evaluate  W(a,l,2jw).  It  should  be  noted, 
however,  that  to  lowest  order  on  both  Aea  and  A,  this  operator  can 
be  evaluated. 


APPENDIX  E 
DIRECT  CCRRELATICN  FINCTION 


In  Appendix  B  we  obtained  the  formal  result, 

K(a,l—«iM)Bjft«,l—«|ft)X*'1(a,l|   )K(a,U  J  (B-10) 

with  the  short  time  limit, 

K(a,l*,,sjof=«|K(a,ljto))=^(a,l',?s;a):^)>C",1(a,lja)=0o)K(a,lsu)). 

The  purpose  of  this  appendix  is  to  examine  the  operator  PC      (a,l;ajs=<») 
which  appears  in  the  operator  V(a,l;u):=«))*  in  the  limit  X    ->0. 
In  this  limit, %  (a,ljur=°°),  defined  in  equation  B-15,  becomes: 

X (a,l|eoS»)L    (l)p~1(a)=  L  [P(a)f  (v  )L    (iJp^aWfdteJpC^) 

ea  a)   .  o     i     e« 


ea*     r  01 


L     (2)p-1(a)-nf  (v  )  d(l)P(a)f  (v  )L    (Dp^a)],  (E-l) 


o     l     ea 


where  the  classical  limit  has  also  been  taken.    By  taking  advantage 
of  the  symmetry  of  L     (l)  we  find. 


ea 


X(a,lja>-)L     (Dp^a):  lp(a)[f(v  )L    (l) 

cd  {0  1        **a 

■*nfd(2)p(a)f  (v2)h(l2)Lea(2)]P"1(a),  (E-2) 

where  we  recall  from  Appendix  D  that 

P(l2)=f  (v   )f  (v   )(l+h(l2)).  (D-ll) 

K  0       10      2 

Equation  E-2  is  not  yet  in  a  convenient  form  to  determine 
its  inverse.     In  order  to  do  this  we  will  Fourier  transform  E-£ 
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FT[X.(a,ljur!»)L(l)p"1(a)  =  Ip(a)[f  (v  )(l*n[dv  h(t)f  (v  )L -  (t))]p"1(a) 

ea  (ij       0   1        z     u   /  <•« 

=  Lp(a)f  (v  )[l+Pnh(S)]L pa(k)p"*1(a)  (E-3) 

a)      0   1     v 


ea 


where  the  operator  P  is  defined  by 


P..F(x,v)= 


dvf  (v)F(r,v). 

0 


(E-4) 


It  is  now  possible  to  solve  E-3  for  L     (£)p~  (a)  with  the  result 


ea 


i-l  -iz-'u-ii 


-li 


LeaU)p*"I(a)=uEl+Pvnh(lt)r1p"1(a)f"1(vi)FT[X(a,lsaF»)Lea(l)p   '(a) 

<=o>(W>    -J31^  p"1(a)f'*1(v  )FT[X(a,ljaF-)L  aU)p"1(a)         (E-6) 
V  l-*nh(lc)  °       x  ea 

where  we  have  used  the  fact  that  P2=P  .     We  now  define  the  direct 

46 
correlation  function      given  by 


nC(Ic): 


Dhjfl 

l-hnh(ic) 


(E-6) 


with  the  inverse  Fourier  transform 


nC(l2)=  t^St  (dfe-^GW  -J^^ 
(2it)     J  l*nh(ic) 


(E-7) 


The  inverse  transform  of  E«*  can  now  be  taken  with  the  result: 


L     (l)p""1(a)=a>(l^)p"1(a)f  -1(v  )*(a,l|aF-)L (Dp'^a) 
ea       H  Oi  ea 


J 


d(2)f  (v  )C(12)F(2). 

0       1 


CF(l)=n 
We  are  now  able  to  rewrite  equation  B-10  in  the  following  form, 


(E-8) 
(E-9) 


K(a,l'"S3ar=«)|K(a,lja!))=^(a,l««»sjaF-)co(l-C)P""1(a)fo(v  )K(a,lj  J. 


(E-10) 


To  be  consistent4%(a#l,,*sjaf=«))  must  be  taken  to  the  same  order 
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-lj 


in  the  coupling  strength  as%      (a,ljaF«)*  hence  E-10  becomes 
K(a,l-»sjaF=o|K(a,l5(o))=fn(vi)fn(vo)[(l+P_)(l+h(l2))-^|d(3)fn(v^) 


0      10      2 
g(l23)P     ](1-C)f"*1(v  )K(a,lSJ. 

'  31  0         1 


21 


0      3 

(E-ll) 


Now  that  we  have  calculated  K(ajl,2jaf=«>|K(a,l5u)))  to  zexoth 
order  in  the  electron  atom  coupling  strength  we  can  now  calculate 
V(a#ljo)==°)  in  the  same  approximation  by  substituting  E-ll  into  B-16 


and  considering  only  the  integral  over  L' 


21 


ee 


r21 


V(a,l;o)=*)=n  d(2)r  'f  (v,  )f  (v  )[(l+P  )(l+h(l2))  + 


ee  o  1  0  2 


21 


+n 


d(3)f  (v  )g(l23)P  ](l-C)f  x(v  ). 

0   3        31       0    1 


(E-12) 


The  term  containing  g(l,2,3)  may  be  rewritten 


.21 


d(2)L^f  (v  )f  (v  )n 
ee  0     !     o     2 


d(3)g(l23)f  (v  )(1-C)f~l(v  )K(a,3;J 

0      3  0        3 


=ftfd(2)[nfd(3)L^f  (v)f  (v  )g(l23)]f  (v  JCl-eJf^Cv  )K(a,2ju)) 

J  J  eeoi03  02  02 

(E-13) 
g(l,2„3)  can  be  eliminated  from  equation  E-13  by  using  the  second 
equation  of  the  equilibrium  hierarchy: 


=«nf  (v  ) 

0      2 


d(3)L^f  (v  )f  (v  )g(l23). 

ee  o     l     o     2 


After  integrating  this  equation  over  v  we  get 


-iv  •$  f  (v  )g(l2Xlf  (v  )g(l2)=Hi  d(3)Lf'f  (v  )f  (v  )g(l23).  (E-14) 


1   1  0   1 


ee  o  l 


31. 


ee  o  l  o  3' 


This  can  now  be  substituted  into  E-13  to  give 


d(2)L    f  (v  )f  (v  )n  d(3)g(l23)f  (v  )(l-C)f"Hv  )K(a,l;J 

ee0i02j  03  03 

Ri[d(2)[-Lee+ivi^i]fo(vi)fo(v2)g(l2)(l-C)fo(v2)k(a,l;to).     (E-15) 
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Equation  E-15  may  now  be  combined  with  E-12  to  yield 
VUjljor^Mnv^f^Vj) 

From  the  definition  of  C,  this  can  be  shown  to  be  equivalent  to 


d(2)f0(v2)g(l2)(l-C)f;1(v2)P2i.       (E-16) 


.-l 


VUljuf^^v^V^v^Cf"1^).  (E-17) 

47 
This  is  the  result  obtained  by  Zwanzig      and  others  for  the  short 

time  kinetic  equation  correcting  the  Vlasov  equation. 


APPENDIX  F 
DYNAMIC  SHIELDING 


The  general  expression  for  the  memory  operator  in  the 
classical  path  approximation  is  given  by 


M(  a))=~in 


d(l)L     (1)  ,-.\  •}      GL     (ljp-^a).  (III-B-28) 

ea         urLu;-V(a,l;J 


In  Chapter  IV  we  observe  that  many  of  the  existing  classical 
approximations  to  the  line  shape  can  be  obtained  by  applying,    some 
weak  coupling  limit  to  this  equation.     Expressions  derived  there 
for  the  denominator  of  equation  III-6-28  have  some  terms  which 
contain  the  permutation  operator,  P. ..     A  convenient  technique  for 
dealing  with  expressions  of  this  type  is  to  use  them  to  form  a 
generalized  dynamically  shielded  electron-atom  interaction.     This 
allows  us  to  rewrite  equation  III-E-28  in  the  form: 


M(o)  )=-in 


,D  /,x i „T     /,x  -l 


d(1)  «(1)  TBtr.  i. )  6  «(1)p  la)       (F-1) 


where  V'(a,l;w)  contains  all  contributions  to  V(a,l;u))  which  are 
not  in  the  form  of  integral  operators,  and  L     (l)  contains  the 
rest. 

From  Section  III-B  we  recall  that  the  memory  operator  can  be 
written,  in  the  classical  path  approximation,  as 


M(w)=n 


d?idVea(*  ^(a,ljj 


K(a,ljJ= ^    ■    ;  - T5Lea(l)P(a)'  (F"2) 

a  \        J  6a 

To  show  that  the  form  for  M(u)  indicated  in  equation  F-l  is  possible, 
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we  will  now  consider  the  operator  V(a,lsw).  We  will  define  an 
operator^  (a,l,2joj)  such  that  V(a,lju))  becomes 


V(a,l;o>)F(a,l;u))^(a,l;u))F(a,l;tJ+  d(2)^(a,l,2;  jF(a,l;w)    (F-3) 

where  F(a,ljw)  is  an  aribtrary  function.  We  will  now  specify  that 
<£(a,l,2jw)  be  a  function  of  x^-x,,  and  v,  only  so  that  we  can 
rewrite  equation  F-2  in  the  form 

V(a,l*))F(a,ljaj)=^r(a,l;u3)F(a,l;o))  + 


dr2^(a,r2-^1,vi,u) 


dv2F(a,r2,v2,w). 


(F-4) 


The  motivation  for  this  restriction,  which  is  valid  for  all  of  the 
approximations  to  V(a,l;w)  considered  in  Chapter  IV,  comes  from 
equation  F-l  in  which  is  is  dv^Ca,^")  for  which  we  must  actually 

■ 

solve.     Then  we  rewrite  equation  IIIH3-28 

fo-L  +iv,^ -L    (l)-V'(a,lsw))K(a,l;co) 
a       i     i    ea 

-  dr2^(a,r2-ri,vljw) 


dv2K(a,r2,v2;w) 


=GL     (l)p"1(a). 
ea 


(F-5) 


We  now  define 


K°(a,r1,v1;aj)= 


1 - GL     (Dp"1  (a) 

w-L  +iv /f-L    (lHr(a,lio>) 


(F-6) 


a       i     i    ea 
which  allows  us  to  rewrite  equation  F-5, 


dVjK^r^Vjjw)-  c 


dv. 


dx. 


or-L  +iv  «V  -L     (l)-\T(a,l;co) 
a       i     i     ea 


^(a^-r^v^w)  dVjKCa,^,^;^)- 


dv^K^r^jio). 


(F-7) 
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For  convenience  in  notation  we  will  define 


Rtar^v  !$<*>)= 


o>-L  +iv,^,-L     (l)-V'(a,l;w) 
a       i     i     ea 


(F-B) 


We  now  Fourier  transform  equation  F-7 


T>    -*■ 


dv jK(a,k,v iiu)*  J7~p 


dv. 


dk1*(a,Jc'-ic,  v,;aj)^(a,Ic^v-|a))x 


dv2K(a,lc^v2sw)=  dv^HaXvjiu). 

This  allows  us  to  define  the  operator 

^(a^lc^v  joj) 


(F-9) 


(F-10) 


so  that  equation  F-9  may  be  written, 


d£"i%aj£»x  so>) 


dv^KUJc^v  jw)- 


dv  K°(a*K*v  |w). 


Now  we  define 

d1c"D""l(aX1c'"sw)D(aX1c';w)=6(M') 


with  the  useful  result 


dVjKCa/^Vjju))^ 


dv^Ca^yv^ju,). 


(F-ll) 


(F-12) 


d1tfr*(aXt'*<i>) 

In  order  to  utilize  equation  F-13  we  note  that  equation  F-l 
can  be  rewritten  in  the  form 


(F-13) 


in 


MU)-  -^p 


d£L    (-IT) 
ea 


d^KCa^^iJ. 


(F-14) 


We  can  now  substitute  F-13  into  F-14  to  get 


M(w)=-  fc   dtc'[  dtcLea(-lc)Er1(a,W'lw)]|dvi«:»(a,k',^;J      (F-15) 
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with  the  inverse  Fourier  transform 

MU)=in[dr  LD  (?  )  [dv  K°(a,r  J  ,w) 


-in 


fd?  I?  (?  ) 


1     l 
i 


.1 


o-L-LjU-L     (lM'-UljJ 


GL     (l)p     (a):   (F-16) 


or 


a    e  ea 


L     (r  ),  the  generalized  dynamically  shielded  electron  atom  interaction, 

63 

is  defined  by 


i?a(*  >>  (^F  ff*^'*"4^  i^C-Mcr^**.* j. 


ea 


(F-17) 


The  form  of  equation  F-15  is  particularly  useful  because  it 

incorporates  all  of  the  effects  of  the  difficult  mean  field  integral 

operators  contained  in   (a,l,2;u))  into  the  dynamic  shielding 

operator*  D  (a,k,k'jio). 

The  shielding  supplied  by  G,  in  the  case  of  the  electron-atom 

interaction  appearing  in  the  numerator,  and  by  V'(a,l;co)  in  the 

denominator  is  easy  to  calculate  in  the  weak  coupling  limit,  and 

is  discussed  in  Chapter  IV.     On  the  other  hand  the  dynamic 

shielding  of  L     (r  )  is  more  difficult  and  requires  a  more 
ea 

complicated  approximation  procedure.     We  start  by  noting  that 
S(a,r  ,v  vjto)  may  be  written; 


F(a,r  .v  :w)= 


l'  l' 


co-L  +iv"  »v        w-L  +iv,  «V, 
ail  all 


— T~i(L    (l)*V(a,lju))x 

63 


xRCa^r^v^co) 

which  in  turn  allows  us  to  rewrite  the  second  term  on  the  left 
hand  side  of  equation  F-*: 


(F-18) 


dv- 


dr    R(a,l  .v  ji,)U(a,i  4  .v  ju) 

2  11  2       11 


dV2K(a,r2,v2;ai): 


dv 


1   w-L  +iv   »V* 
a       i     i 


^-((a^-ff^Vjju) 


dv2K(a,r2,v2Jto)-i 


dv. 


dr 


1   u-L  +Iv   *$ 
a       l     i 


(Loa(*1)'*V'(a,ls««j))R(a,r1-vju)*5C(a,x,  •«.!«) 


This  can  be  Fourier  transformed  with  the  result 
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dv2K(a,r2v2jw).  (F-19) 


FT 


d? 


dx  RCa^r^v^uJ^Ca^-tf^v^u)  dv2K(a,r2,v2;w) 


*T    '   gC(a»k,v  ;a>)  dv,K(a*k*v-|u)- 75*rr 

+^    .£   ^  1  J         2  2  V2!"         J 


dv 


ur*L 


dvxdk  dk 


a     l 


1 L    (le-t")R(a»ic"-k',v1ia»)i(a,k>1|w) 

u>-L  +v>Ic     ea  1 


dv2K(a,ic',v2sw), 
(F-20) 


Coirbining  this  result  with  equation  F-8,  F-9,   and  F-10  suggests 
that  D(a,k,k'jto)  should  be  divided  into  two  parts:     one  which  is 
local  in  k«space  and  one  which  is  nonlocal* 


:(o) 


D(a,lU'l  0))=OVU;(a,kMc'j  ud*Du'(aXk 'i  ti 


(U, 


D(0)(aX^J=6(^')ll+i 


dv,  i~— <£Uk~;u))]=6(k-k')e(kju)-L  ) 

^  -L  ^  .£  a 

a     l 


D(l)(a,lU'jJ=  T^Ts-fdv  dlt *— -  L    («')R(a,t"-Jc',v1|«) 

l27T;  J      l  or-L  +v  4    ea 

a     l 

x^(a,^vi5J.  (F-21) 


If  we  assume  that  e(£jar-L  )  the  mean  field  part  of  V(a,l;to)  is  the 

3 

Vlasov  operator,  then  e(kjar«L  )  is  the  usual  frequency  dependent 


dielectric  constant.  .  We  now  have: 


■n(0) 


(1), 


dlTirHa,ir,ir-sJlDvu/(a*^^jw)+DVi/(a,P%ic'jU)]=6(ic-P)   (F-22) 


or,  using  equation  F-£l 


D"1  ( *,tMt  %  J  e(kj  uj-L  )  +  fdlt^D"1  ( a,£,)t'  'i  to) 


D(l)(a,lt->lt'jJ=6(ir'-kv) 


(F-23) 
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ff-Cfct-id-  -dp*1  -    r*     ;   x 

e(kjor-L  )       e(kjorL  ) 
'«fdc'-6"1(a.M"l«)B(1,(«.t".*'ltt).  C*"*4) 

Thus,  using  equation  F-16,  we  have  a  formal  expression  fox  the 
dynamically  shielded  electxon-atom  interaction.     It  should  be 
observed  that  D^     (a*k",,k';io)  does  not  contribute  to  lowest  order 
in  the  electron-atom  coupling  strength  and  equation  F-16  becomes 

LD  tfj*  tAt  ke"lk'r    -« •  (F-25) 

3 
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